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Nomenclature 


[A] apparent mass matrix, see Eq. (90) 

[B] damping matrix, see subscripts below for definitions of components 

[B] dimensionless stiffness matrix, see Eq. (89) for component scale factors 

C clearance (cylindrical seal) or reference film thickness (face seal) 

e x ,e z displacements of center of seal in x,z directions, see Fig. 1 
e x ,e z dimensionless displacements, e x /C, e z /C 

e x0 ,e z0 dimensionless displacement amplitudes, see Eq. (60) 

F function defining derivatives, see Eqs. (56)-(58) 

f shear factor (% friction factor) 
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l f face seal flag parameter, l f =1 for face seal, l,=0 for cylindrical seal 
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1 , 7 unit vectors in 0,s directions 

S(Z) imaginary part of complex number, Z 

J k see Eq. (64) 

[K] stiffness matrix, see subscripts below for definitions of components 

[K] dimensionless stiffness matrix, see Eq. (88) for component definitions 
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derivative matrix used in numerical solution, see Eq. (58) for definitions of components 

seal length, see Fig. 1 

applied moments about x, y axes 

dimensionless moments, M x /(p 0 rjj), ty/l /(p 3 f ) 
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pressure 

dimensionless pressure, p/p 0 
second order pressure 
dimensionless pressure, p'/p 0 
reference pressure 

dimensionless parameter, pV 0 r 0 /(4C 2 p 0 ) 

pressure at upstream side of jump in film thickness 

second order pressure at upstream side of jump in film thickness 

component of second order pressure, p', associated with 5,,, see Eq. (65) 

complex, second order, dimensionless component pressures, see Eqs. (66)-(67) 

complex S dependent factors of pS £nd p respectively, see Eq. (67) 

pressure at exit side of seal 

pressure at inlet side of seal 

dimensionless pressures, p ex /p 0 , p jn /p 0 

local pressure gradients over grooves in tangential (p g 9 ) and transverse (p g s ) directions 
dimensionless local pressure gradients r 0 p g6 /p 0 , r 0 p gs /pt>, see Eqs. (24)-(25) 
local pressure gradients over ridges in tangential (p r8 ) and transverse (p rs ) directions 
dimensionless local pressure gradients r 0 p r6 /p 0 , r 0 p rs /p 0 , see Eqs. (26)-(27) 
power loss 

dimensionless power loss, P/(Cp r 0 <§) 
volumetric flow rate in transverse (s) direction 

viii 
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q g9 ,q gs dimensionless flow components over grooves, u g h g , v g hg 

q^,q re dimensionless flow components over ridges, u r h r , v r h r 

q n flow rate in n p direction relative to spiral groove motion 

q n dimensionless flow rate normal to spiral grooves, qn/(hoV 0 ) 

R Reynolds number, see Eq. (9) 

R* characteristic inertia parameter, 2CR/r 0 

R a local Reynolds number associated with surface (a), see Eq. (6) 

R b local Reynolds number associated with surface (b), see Eq. (6) 

r radial coordinate, taken as r 0 for cylindrical seal 

r dimensionless coordinate, r/r „ taken as 1 for cylindrical seal 

r 0 reference radius, taken as outside radius for face seal and shaft radius for cylindrical seal 

3t(Z) real part of complex number, Z 

S dimensionless coordinate, s/r 0 

Sj dimensionless coordinate s/r 0 

S L ,S R dimensionless boundary coordinates s L /r 0 , Sr / r 0 

S in ,S ex dimensionless coordinates s in /r 0 , s ex /r 0 

s transverse coordinate, s = r for a face seal and s = z for a cylindrical seal 
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s L left boundary coordinate, inside radius for face seal and s,_ = -L/2 for cylindrical seal 

s R right boundary coordinate, Sr = r 0 for face seal and s R = U2 for cylindrical seal 

s ex s coordinate at exit from seal 

s in s coordinate at inlet to seal 

t time 

t dimensionless time, V Vr 0 

~t p unit vector tangent to grooves, see Fig. 4 

U first order tangential (0) velocity component 

0 dimensionless velocity, U/V 0 

u tangential (0) velocity component 

u dimensionless velocity, u/V 0 

u' second order tangential (0) velocity component 

u' dimensionless velocity u'/V 0 

u fluid velocity vector 

u a velocity of surface a, see Fig. 2 

u b velocity of surface b, see Fig. 2 
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w x ,w y ,w z 

W x ,w v ,w z 


^ynewj 


Ap 

Ap 

Ap' 

APg,Ap r 

APg,Ap r 


dimensionless velocity, u in /V 0 

component of u', associated with 5 k , see Eq. (65) 

complex, second order, dimensionless component tangential velocities, see Eqs. (66)-(67) 

complex S dependent factors of G ^nd u respectively, see Eq. (67) 

first order transverse (s) velocity component 

dimensionless velocity, VA/ 0 

reference velocity 

transverse (s) velocity component 

dimensionless velocity component, vA/ 0 

second order transverse (s) velocity component 

dimensionless velocity v'A/ 0 

transverse velocity over grooves (v g ), ridges (v r ) 

dimensionless velocity components v g A/ 0 , v r A/ 0 

s velocity component at upstream side of jump in film thickness 

dimensionless velocity VjA/ 0 

second order s velocity component at upstream side of jump in film thickness 
dimensionless velocity 

component of v', associated with 5 k , see Eq. (65) 

complex, second order, dimensionless component transverse velocities, see Eqs. (66)-(67) 

complex S dependent factors of $ >and v Respectively, see Eq. (67) 

first order axial load applied to face seal 

dimensionless load, W z /(p 0 r^) 

applied loads in x, y, z directions 

dimensionless loads, w k /(p 0 i^), yv /(p^r ), vy /(g f \ 

coordinate variables, see Fig. 1 

column vector of independent variables at old grid point, see Eqs. (56)-(58) 
column vector of independent variables at new grid point, see Eqs. (56)-(58) 
groove to pitch ratio, A0 g /(A0 g +A0 r ) 

spiral groove angle, angle between grooves and surface velocity 

jump in film thickness (downstream - upstream) 

dimensionless jump in film thickness, Ah/C 

jump in pressure at Sj, p-pj 

dimensionless jump in pressure, Ap/p 0 

dimensionless second order jump in pressure at Sj, (p'-pj)/p 0 

spiral grv. downstream - upstream pressure jump at entrance to grooves, (Ap g ) or ridges, (Ap r ) 
dimensionless spiral groove pressure jumps, Ap g /p 0 , Ap r /p 0 

variable increment between grid points used in numerical solutions, see Eqs. (57)-(58) 
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A0 circumferential extent of groove-ridge pair, see Fig. 4 

A9 g circumferential extent of groove, see Fig. 4 

A0 r circumferential extent of ridge, see Fig. 4 

5 groove depth, h g -h r 

5 dimensionless groove depth, 5/C 

5 k dimensionless component film displacements, see Eqs. (59)-(61 ) 

5 k ,?jt complex dimensionless component displacements associated with + Q and -ft resp., see Eq. (63) 

e k set of amplitudes, see Eq. (64) 

t, loss coefficient due to contraction 

H k see Eq. (64) 

0 angular coordinate, see Fig. 1 

A correction factor applied to t c , see Eq. (94) 

p viscosity 

J; loss coefficient, see Eq. (1 8) 

p density 

t tangential shear stress 

t dimensionless shear stress, r j/(Cp \ 

T a tangential shear stress, surface (a) 

T a dimensionless shear stress, r 0 T a /(Cp 0 ) 

T a shear traction vector for surface (a), see Fig. 2 

T b tangential shear stress, surface (b) 

T b dimensionless shear stress, r 0 T b /(Cp 0 ) 

T b shear traction vector for surface (b), see Fig. 2 

t c dimensionless tangential Couette shear stress, (f a + f b )/2 

t c Couette shear traction vector, ( T a + t„)/2 

t p dimensionless tangential Poiseuille shear stress, (f a - i b )/2 

t p Poiseuille shear traction vector, (T a - f b )/2 

t dimensionless effective shear stress, see Eq. (92) 

T* g f values of t obtained from first order solutions for grooves and ridges respectively 

turbulent shear functions, see Eqs. (1 4)-(1 5) 

(tf.ty* global functions defined by Eqs. (42)-(43) 

tJ>g,*4J g turbulent shear functions over grooves, see Eq. (22) 

<t> r ,4J r turbulent shear functions over ridges, see Eq. (23) 

X function for predicting pressure change at jump in film thickness, see Eq. (17) 

4 J rotation about y axis, see Fig. 1 

qj 0 dimensionless amplitude, see Eq. (60) 

Q angular velocity of oscillation 
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0 dimensionless angular velocity, Qr /y 0 

0 + ,Q- seeEq. (71) 

to angular velocity of rotor 

a) dimensionless angular velocity, cor /y 0 

Subscripts 

a,b relating to surface (a),(b) shown in Fig. 2, generally (a) is moving, (b) is stationary 

ex exit 

g grooves 

H,U,V partial derivatives with respect to H, 0 and V respectively 
in inlet 

J jump discontinuity in global film thickness 

k index associated with displacements, see Eq. (64) 

L left or lower value of s 

R right or upper value of s 

r ridges 

x pertains to force or displacement in the x direction 

y pertains to force or displacement in the y direction 

4> pertains to moment or rotation about x axis 

HJ pertains to moment or rotation about y axis 
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1. Introduction 


Spiral groove bearings and seals are used to provide stability, load support and pumping for both cylindrical 
and face seal geometries. In the case of a cylindrical bearing, grooves are usually designed to pump against 
each other in a symmetric arrangement to provide enhanced stability. A lightly loaded cylindrical seal 
operating at a low axial flow rate will produce a force that is nearly 90 degrees out of phase with the 
displacement which will tend to destabilize the rotating shaft. The introduction of spiral grooves can 
significantly increase the component of force in phase with the displacement and decrease the out of phase 
component thereby improving stability. 

In the case of a face seal or thrust bearing, spiral grooves are often introduced as the primary means of load 
support. Since a symmetric arrangement is not possible in a radial geometry, the grooves are usually 
designed to pump towards an ungrooved dam region. The resistance of the dam region increases as the film 
thickness decreases hence the pumping pressure rise increases thereby giving rise to a positive axial 
stiffness. The spiral grooves can also be used to pump against an applied pressure gradient thereby resulting 
in either reduced or reversed leakage. 

The computer codes developed here are oriented toward the prediction of performance characteristics of 
incompressible cylindrical and face seals with or without the inclusion of spiral grooves. Performance 
characteristics include load capacity (for face seals), leakage flow, power requirements and dynamic 
characteristics in the form of stiffness, damping and apparent mass coefficients in 4 degrees of freedom for 
cylindrical seals and 3 degrees of freedom for face seals. These performance characteristics are computed 
as functions of seal and groove geometry, load or film thickness, running and disturbance speeds, fluid 
viscosity, and boundary pressures. 

The basic assumptions that have gone into the computer code are listed below: 

1 . The flow is assumed to be isothermal and incompressible. 

2. Turbulence is treated with an extended form of the Hirs bulk flow model (reference 1), generalized 
to include separate and arbitrary friction factor - Reynolds number relationships for each surface. 

3. Inertia effects (which throughout this report will refer to additional effects to those inherent in 
turbulence) associated with film discontinuities will be treated with the use of loss coefficients. 

4. Circumferential and transient effects are treated with the use of small perturbations to a steady state 
first order solution for a concentric, aligned seal. These effects are characterized by either by 
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stiffness and damping coefficients that are dependent on the disturbance frequencies or by stiffness, 
damping, and apparent mass coefficients. 

5. The film thickness is assumed to be small compared with seal lengths and diameters but large 
compared with surface roughness. 

.6. Narrow groove theory is used to characterize the effects of grooves by a global pressure distribution 
without requiring computations on a groove by groove basis. This has previously involved neglecting 
edge effects and local inertia effects associated with groove to groove pressure variations. In 
general, narrow groove theory is valid when there are a sufficiently large number of grooves so that 
2n:/N g « 1, where N g is the number of grooves. This theory has been extended to account for the 
pressure jumps arising from the effects of local inertia associated with the sudden changes in cross 
section at the inlet to each groove and ridge which can be of the order of N g larger than the other local 
inertia effects which have been neglected. 

7. Although transverse (axial for cylindrical seals or radial for face seals) variations in the film thickness 
profiles are permitted, machined surfaces for seals are assumed to be axisymmetric with the 
exception of effects of spiral grooves. 

The above assumptions still leave the code applicable to a broad range of cylindrical and face seal 
applications. Practical designs should contain a fairly large number of grooves to ensure smooth, isotropic 
operation. Elastic and thermal distortions as well as machining tolerances should also be estimated to validate 
the axisymmetric clearance assumption and may be characterized in part by inclusion of film variations in the 
transverse direction. The overall accuracy of the program will depend on the grid size used. Factors such 
as large film thickness variations in the transverse direction, inclusion of very small transverse inertia effects 
and large values of the length to diameter ratio could require an increased number of grid points. 

A derivation of the equations governing the performance of turbulent, incompressible, spiral groove cylindrical 
and face seals along with a description of their solution is given in the next section. This will be followed by 
a description of the computer codes including an input description, sample cases and comparisons with results 
of other codes. 
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2. Formulation and Method of Solution 


Coordinate variables will be used to make the equations developed here applicable to both cylindrical and face 
seals and may be defined with the aid of Figure 1. The circumferential coordinate, 0, is as shown in Figure 
1 . The transverse coordinate is described by the variable, s, which is taken to equal the radial coordinate, 
r, for a face seal and the axial coordinate, z, for a cylindrical seal. The quantity r, when it appears will denote 
radial position for a face seal and should be set equal to the shaft radius, R<, for a cylindrical seal. The 
symbols Tj will be used to denote unit vectors in the 0 and s directions respectively. 


Formulation of bulk flow equations for turbulent lubrication 


The equations for turbulent bulk flow to be presented here will be based on a generalization of the Hirs model. 
Figure 2 shows the velocities and forces associated with a differential element in the 0 direction. The bulk flow 
integrated momentum equations for the 0 and s directions are 


ph 


ph 
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du 

4- \/ 

u 

-i- 

3u uvl 

+ 1 

_ h 

3p 

3t 
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ds 
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30 r J 
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+ v — 

u 

+ — 

3v _ u 2 l f 

= -h 

3p 
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ds 

r 

30 r j 


ds 


b a ‘ 


“b a' 


The integrated continuity equation is 


13... 13,,. 3h 

(rvh) + (uh) + — =0 

r 3s r 30 3t 


(1) 

(2) 


(3) 


In the above equations u and v represent the components of the velocity vector, u = u i + v j , p denotes the 
pressure, x a and x b are the shear stresses acting on the fluid at surfaces (a) and (b) respectively as shown 
in Figure 2 and l f is a flag parameter set equal to 1 for a face seal and 0 for a cylindrical seal. 


Each of these shear stresses is assumed to act in the direction of the fluid velocity relative to the respective 
surface and is related to that velocity by means of a coefficient of resistance or friction factor as follows: 


P I u - u a | f a 


2hp | u - u. 


(u -u a ) 


4 h 


u.) 


(4) 
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Figure 1 Coordinate system for spiral groove analysis. 
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Figure 2 Velocities and forces on a differential element in the 0 direction, 

for surface (a) and for surface (b) 


T b = --Pl u 


u J f b 


2hp | u - u t 


(5-u„> = -7^R b f b (R»)(u -u b ) 

4 h 


(5) 


where the shear factors f a and f b are one fourth of the resistance coefficients or friction factors for the 
respective surfaces (reference 2) and are functions of the respective Reynolds numbers R a and R b which are 
in turn defined as 


R a = 2h | u - u a | p/p , R b = 2h | u - u b | p/p . 


(6) 


The use of various shear functions in the prediction of turbulent behavior of bearings and seals may be found 
in references 3-6. Power law relationships of the form f = nR m have been used in references 3 - 5 to 
characterize both smooth and rough surfaces. If both surfaces are smooth the relationship 


f.(R.) = n o R £ 


n o^b 


(7) 
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may be used with n 0 and m 0 obtained from the Biasius Equation (reference 1) as n 0 = 0.0751, m 0 = -.25 for 
turbulent flow (Reynolds numbers ranging from 4,000 to 100,000). The corresponding values for laminar flow 
are n 0 = 24 and m 0 = -1 


Equations (1) - (3) may be expressed in dimensionless form with the use of the following dimensionless 
variables 


h 


h_ 

i 

C 



S = 




x 


liJL 

C Po 


(8) 


and parameters 

M 




P* 


P V o r o 

4C 2 p 0 


( 9 ) 


In the above equations V„ represents a characteristic velocity, p 0 a characteristic pressure, r 0 and C represent 
the radius and clearance for a cylindrical seal and the outside radius and characteristic film thickness for a 
face seal. 

Henceforth surface (b) will be taken to be stationary and surface (a) will be taken to move with a 
circumferential velocity rw. The Reynolds numbers R a and R b given by Equation (6) then become 

R a =Rh \/(0 -rco) 2 + v 2 , R b =Rh v / u 2 +v 2 . (10) 


The integrated momentum and continuity equations may now be expressed in the dimensionless form 


1 dp 
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( 13 ) 


1 d(rvh) + 1 d(uh) + dh 
f dS r 30 pH 


where 


a>(u,v,h) 


(Q-f«)RA(R a ) +QR b f b( R b) 


and 


Y(u,v,h) 


Rafa(Ra) + Rby R b) , 


(14) 


(15) 


Equations (11) - (13) together with the functions defined by Equations (10), (14) and (15) constitute the 
dimensionless differential equations for the generalized Hirs bulk flow model under consideration. These 
equations require a prescribed film thickness profile, h(S,0,t), definitions of f a and f b such as those given by 
Equation (7) and boundary conditions which will be treated below. 

Boundary and continuity conditions 

This analysis is limited to essentially axisymmetric surface shapes and constant boundary pressures. 
Although spiral grooves will be considered, the groove angle, depth and spacing will be independent of 0. 
Individual pads and partial films are not treated here. As such, all 0 (and time) dependence will be associated 
with displacements of the sealing surfaces. The only conditions on 0 are the requirement of periodicity. That 
is, pressure and flow rates at 0 + 2n must equal those at 0 for all 0. 

Variations of surface shapes in the s direction, including discontinuities are considered. The upstream values 
of v and p at discontinuities in h are denoted by Vj and pj as shown below. The film thickness changes from 
the upstream value of h - Ah to the downstream value of h across the jump. 


p , v , h - Ah 
' j 

p, v, h 

-Ah 

s flow direction 


Figure 3 Jump in film thickness. 
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Since the flow rate normal to the jump must be continuous 


(h -Ah)Vj = hv @ s = Sj , 
where Sj denotes the s coordinate at the jump. 

The pressure change at Sj may be expressed in terms of a loss coefficient based on downstream conditions 
(reference 7) as follows: 


Pj + -jP v j 2 = P + ^pv 2 (1 +0 @ s = Sj 


The preceding two equations may be rearranged and expressed in dimensionless form as 


v = (1 )Vj , Ap = p - pj = lp*R*x(v,h,Ah) @S = Sj , 

h 2 


( 16 ) 


where 


X(v,h,Ah) = 




2 

1 - 

h 

+ 5 

, h - Ah , 




v 2 


( 17 ) 


The loss coefficient ij, due to a contraction will in general be a function of the axial Reynolds number and will 
be denoted by the function £(Rhv) but is frequently taken to be constant (references 3 - 5). The loss 
coefficient for an expansion may be found in reference 7 and is given together with the corresponding 
relationship for a contraction below: 


C(R,h,v) 



1 


h 


\ 2 


h - Ah ) 


Ah < 0 (contraction) 
Ah > 0 (expansion) 


( 18 ) 


The above considerations together with the requirement that the tangential component of velocity is 
continuous at jumps across constant s boundaries provides sufficient information to prescribe the boundary 
conditions at the inlet and exit to the seal as well as the continuity conditions at jumps. 


The tangential component of the velocity, u in and the pressure, p in are prescribed at the inlet side of the seal. 
The relationship between p in and v at the inlet is obtained by setting p d = p in and evaluating % as Ah - -«> in 
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Equations (16) - (18). These conditions may be written in dimensionless form as 


u ; . 


Pm + -jP*R*X(v,h.-°°) = p in 


1 s* 


p*R* (1 + C)v 


r,2 


s =s : . 


(19) 


The pressure at the exit from the seal will be taken to be continuous and equal to the exit side pressure, p ex . 

P = Pex @ S = S ex ■ (20) 

It should be noted that S in may correspond to either the left or right hand sides of the seal (lower or upper 
values of S) depending on the direction of flow. S ex will be the value of S at the opposite side of the seal from 

s = s in . 

Narrow groove theory for spiral groove seals 

The development of the narrow groove equations to be presented here will expand upon equations used in 
reference 8 which, in turn, contains a reformulation of the lubrication equations for laminar, compressible 
spiral groove seals presented in reference 9. 

The purpose of narrow groove theory is to obtain equations and solutions for continuous, smooth global 
pressure and velocity distributions that approximate the limiting form of the solution to Equations (10) - (15) 
as the number of grooves becomes large, with the groove angle, p and the groove to pitch ratio, a, held 
constant. This process greatly reduces the complexity of the solution and eliminates the need for the large 
grid that would be required to adequately describe all of the grooves individually. 

The discontinuities in film thickness associated with the grooves will give rise to discontinuities in the 
pressures, pressure gradients, and velocities at the ridge-groove interfaces. The relative discontinuities in 
the velocity components normal to the grooves will be of the same order as the relative film discontinuities. 
The predicted discontinuities in the local pressures at the interfaces would be of the order of the inertia (p‘R*) 
terms on the right hand sides of Equations (11) and (12). It is important to note however, that as both the 
number of grooves and the inertia terms become large, flows will not be able to fully develop within or between 
grooves. 

Figure 4 shows three conceptual states of flow development in the order of increasing significance of local 
fluid inertia and decreasing degree of flow development. These are Case a, where the local pressure changes 
at the inlet to each groove (or ridge) are negligible; Case b, where a wake forms at the inlet to a groove whose 
extent is small compared to the length of a groove and Case c, where a jet forms and prevails over most or 
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all of the length of the groove. These three cases are discussed individually below: 



a. Negligible inlet pressure changes. Spiral groove seals frequently operate at relatively small 
clearances (as opposed to labyrinth seals or even damping seals) and will pump at relatively low flow 
rates compared with those produced by imposed pressure gradients of sufficient magnitude to 
produce large inertia effects. Under these conditions, the flow can be treated as fully developed over 
the entire groove and inertia effects may be neglected in the analysis of the local pressures and 
velocities. All inertia effects in Equations (11) - (20) may, however, be included in the treatment of 
the global pressure and velocities. As such, pumping affected by turbulent shear in regions 
containing spiral grooves would be accounted for and global inertia effects due to circular steps and 
grooves and other axial clearance variations would still be included in all regions. 
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Figure 5 Schematic of spiral groove parameters, global and local pressures, Case a 

A schematic showing a global and local pressure distribution for this case along with some of the 
spiral groove film variables under consideration is presented in Figure 5. The local pressure profile 
is shown by the sawtooth lines whose lower vertices, for purposes of illustration, are connected by 
the global pressure profile, p. The global pressure profile does not necessarily lie at the lower 
vertices of the local pressure profile but could lie anywhere between the lower and upper vertices. 
In the limit as the number of grooves becomes large, curves connecting the upper and lower vertices 
will approach each other (subject to the approximations discussed above). This limiting behavior is 
not true of u, v, h or the local pressure gradients which will have different values over lands and 
grooves no matter how large the number of grooves. 

b. Fully developed films with significant inlet pressure changes. This case extends Case a by 
including local pressure changes at the inlet to each groove and ridge with the use of loss coefficients 
in a manner analogous to that used at the inlets to each region given by Equations (16) - (18). The 
treatment of these pressure changes as boundary pressure discontinuities assume that the length 
of the wake shown in Figure 4b is small compared with other relevant circumferential and transverse 
length dimensions. The relevant length for comparison with the wake length at the start of a region 
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would be the transverse length of the region. In the case of spiral grooves within a region, the 
relevant lengths would be the widths of the grooves and ridges. 



The local pressures shown in Figure 6 are a generalization of those in Figure 5 to allow for local 
pressure jumps at the inlet to each groove and ridge. Unlike Case a, the upper and lower envelopes 
of the local pressure curves will not coalesce as the number of grooves becomes large since the 
magnitudes of the groove and ridge widths do not enter directly into the jump relationships given by 
Equations (16) - (18) and the magnitudes of the jumps will not approach 0 as long as Case b 
conditions prevail. The local pressure jumps for this case will be shown to give rise to a term of order 
(p'R'Ng) which will be retained in the analysis. As with case 1, the other inertia terms will be 
neglected in the local analysis but maintained in the global analysis. These assumptions will thus give 
rise to asymptotic relationships as the number of grooves becomes very large in contrast to the 
limiting relationships corresponding to the assumptions used in Case a. 

c. Jet extends over entire groove. As the number of grooves becomes very large a limiting flow 
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should in fact occur as a jet forms that extends across the entire groove. The bulk flow theory which 
is based on the assumptions that the shear forces on the turbulent core arise directly from interactions 
between the core and the bounding surfaces break down at this point. A rigorous analysis would 
involve consideration of the jet, the wake, interactions between jet and wake, interactions between 
jet and wall etc. and would require three dimensional solutions to the momentum and continuity 
equations which is beyond the scope of the present effort. An alternative approach would involve the 
use of additional empirical relationships to characterize these interactions so that limiting behavior can 
be obtained. Examples of the use of such relationships are given in reference 10 for labyrinth seals 
and in reference 1 1 for straight through screw seals. 

The narrow groove analysis contained in the rest of this section will refer to Case b unless otherwise specified. 
The theory requires the development of expressions for the local quantities in terms of continuous global 
quantities that become valid as the number of grooves becomes large. Local quantities will be denoted by 
the subscript g when referring to grooves and r when referring to ridges. Variables such as p,u,v appearing 
without r or g subscripts will henceforth denote global quantities. The subscript r denoting ridges has been 
used for consistency with existing literature and should not be confused with the radial position variable r 
which will not be used as a subscript. 


Spiral groove patterns for a given region are characterized by three parameters which are defined in terms 
of quantities shown in Figure 6. These are the groove depth 8 = h g - h r , the groove to pitch ratio 
a = A0g/(A0 g +A0 r ) and the groove angle p defined as the angle between the groove and the surface velocity 
direction, measured from the groove. 


One may now introduce the dimensionless groove and ridge flows in the 0 and s directions 


q g e = G- h 


g g 


q„ e = v h 

^gs g g 


q r e = > 3r 


v r h r 


( 21 ) 


and express the turbulent shear functions $ and *F defined by Equations (14) and (15) for the grooves as 


$ = ) , y = Y(^,^i,h g ) 


( 22 ) 


and for the ridges as 


<D r = <D(^>,h r ) , Y r = Y(^>,h r ) . 

h r h r h r h r 


(23) 
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One may now write Equations (11) and (12) for the groove regions as 


Pg,0 


_ 1 dp 

7 30 


" P9,s= ”as 


and for the ridge regions as 


1 dp, 

■ Pr ' e= ~7ie' 


-t 

r ' s 3s' 


p*<I> g + p*R* 


+ v — + 




p*T g + p*R’ 


p*$ r + p*R* 


+ v — + 


p*Y r + p*R* 


+ v 


u 3u 

-|- 

,«) 

r 30 

r 


/ 

u 3v 


7 30 

■f . 

r 

> 

u 3u + 

__ \ 
i ^ 

730 

■f . 

r J 

i 3v _ - 



3S r 30 f r 


( 24 ) 


(25) 


(26) 


(27) 


The unit tangent and normal vectors for the logarithmic spiral as shown in Figure 6 are given by 

t p = cos Pi +sinPj , n p = sin Pi — sin p j (28) 

If local inertia effects were neglected, the pressure would be continuous at each groove-ridge interface and 
the derivative of the pressure in the direction tangential to the interface, Vp*? p , would also be continuous. 
In keeping with the assumption that the only local inertia effect to be considered is the pressure jump normal 
to the interface, continuity of the tangential pressure derivative will be maintained even when this effect is 
considered. Continuity of the tangential derivative of pressure at the interface may be expressed as 

cosp p g 0 + sinp p g s = cosp p r 0 + sinp p rs . (29) 

One may now substitute Equations (24) - (27) for the corresponding pressure derivatives in the above 
equation to obtain 

cosP<& g + sinPTg = cosp <E> r + sinPT r . (30) 
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The flow rate relative to a groove-ridge interface in the direction normal to that interface, q„, must be 
continuous. For grooves on the stationary surface (b), that flow rate would be hu«n p and for grooves on the 
moving surface (a), it would be h(u - rco i )«n p . One may use the dimensionless groove and ridge flow rates 
defined by Equation (21) and perform the indicated dot products to express this continuity condition as 

% = (q g e -fwh g IJsinp - q gs cosp = (q r0 - ffi>h r IJsin p - q rs cosp , (31) 


where l u is equal to 1 if the grooves are on the moving surface and 0 if the grooves are on the stationary 
surface. The two equalities on the right may be rearranged to obtain the relationship 


sin Pq g e - cospq gs - sinpq r9 + cospq rs = fwSl^sinp . 


(32) 


One may now equate flow the across a groove-ridge pair in the 0 direction to the global flow in that direction 
and express the result in dimensionless form as 


a q g e + (1 -°0q r e = • 


(33) 


A similar balance in the s direction results in the relationship 

a q gs + (i -°0q r s = vh . 


(34) 


The global film thickness, h is the value which results in the same flow area as that produced by the local film 
thickness variation and may be expressed in dimensionless form as 

h = ah g + (1 - a)h r = h r + a& (35) 

Equations (30) and (32) - (34) represent four algebraic equations which may be solved for the four local flow 
rates q g9 , q gs , q rt and Iterative solutions are obtained numerically with the use of a Newton-Raphson 
procedure to linearize Equation (30). 

The pressure changes (downstream - upstream) incurred at the entrance to a groove, Ap g , and at the 
entrance to a ridge, Ap r , may be expressed in dimensionless form with the aid of the jump function %, defined 
by Equations (16) - (18), as follows: 
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A P C 


On r 


P‘R*X(-.h g .6) . Ap r = lp*R* X (^,h r ,-8) 


(36) 


The normal flow component, q n , may be expressed in terms of the global velocities and film thickness with 
the use of Equations (31) and (33) - (35) as 

q n = (Oh - rcohl^sin (3 - vhcosP . (37) 


By equating the change in pressure over a groove-ridge pair to the changes in the global pressure, first in the 
0 direction, as shown in Figure 6, and then in the s direction and expressing the result in dimensionless form, 
the following equations are obtained: 


1 dp 

7 30 


“Pg.e + (1 -a)P r ,e 


A Pg + Ap r . - - 

+ fAe ^Kv) 


(38) 


— = ap + (1 -a)p + — ^ — sign(q ryj ) 

3S 9 ' s r ' s FA0 1 tan P| " p 


(39) 


where the sign functions appearing in the above equations are either 1 or -1 depending on the direction of the 
normal component of flow with respect to the corresponding coordinate direction. 


Finally, noting that A0 = 27t/N g and the definition of n p given by Equation (28), one may substitute Equations 
(24) - (27) for the corresponding local pressure gradients and Equation (34) for the pressure jumps in 
Equations (38) and (39), and write the resulting equations in the form 


1 dp 

7 30 


p*^(G,v,h,l f ) + p*R* 


' 3u _ 30 u 30 ^ 

— + v — + 

3S r 30 


3t 


(40) 


3p 

3S 


p*T*(0 t v,h,l f ) + p*R* 


1 3v _ 3v 0 3v ^ 
— + v — + 

3S r 30 


at 


where 


uv 


® = a3> g + (1 -a)3> r + R * 1l f — 


47ir 


rsign(q„sinp) 


X(^.h g ,6) + x (^.h r .-6) 


(41) 


(42) 
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and 


T* = oT + (1 - a)Y r + R* 


-. f °i + — ^ 

r 4tcF | tan P | 


g 


sign(q n cosp) 


X(T L .h g ,5) + x(^-,h r , -8) 

h„ y h r 


( 43 ) 


The Coriolis and centrifugal inertia terms for face seals and the spiral groove local inertia term have been 
combined with the turbulent shear resistance terms in Equations (42) and (43). This has been done for 
convenience because of their algebraic nature and will simplify the development and solution of the first and 
second order flow equations to be discussed in the next section. 


The actual number of grooves, N g , does not usually appear in a narrow groove analysis. It enters here only 
to characterize the cumulative effects of losses due to sudden changes in cross section between ridges and 
grooves. Formally setting N g =0 would therefore only eliminate the inclusion of these local inertia effects and 
would maintain all of the other spiral groove effects under consideration. Caution should be exercised when 
groove or ridge widths become the same order of magnitude as the clearance. This would invalidate one of 
the basic geometric assumptions of the bulk flow theory used throughout this analysis and could have a 
particularly pronounced affect on the validity of the above local inertia terms. 

Equations (40) and (41) represent the generalization of Equations (11) and (12) to spiral groove seals and 
reduce to Equations (11) and (12) when both a and N g are set equal to 0.. The continuity equation 
[Equation(13)], and the continuity and boundary conditions given by Equations (16) - (20) remain unchanged 
subject to the interpretation of u, v, p and h as their global values when a > 0. 

Development of first and second order equations for small displacements 

The first order variables will be the global film thickness, pressure and tangential and transverse velocities for 
an undisplaced shaft. They will depend only on s and their dimensionless forms will be denoted by H(S), P(S), 
U(S) and V(S) respectively. Substitution of the above variables for h, p, u,.and v respectively, in Equations 
(40), (41), and (13) yields the following set of ordinary differential equations: 


<3? (U,V,H,l f ) + R* V 


dU 

dS 


= 0 , 


( 44 ) 


= p‘Y*(U,V,H,l f ) + p* R* V — — 
dS dS 


( 45 ) 


dfrHV) w Q 
dS 


( 46 ) 
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l CL 


The boundary and jump conditions for the first order variables are prescribed by Equations (16) - (20) with h, 
, u, and v replaced by the corresponding first order variables. The boundary conditions are given by 


U = Q,. 


P = Pin + jP*R*X(V,H,-oo) 


Pin - jP’R’O + C)V @s=s in (47) 


and 


P 


S =S, 


(48) 


One may now add a small displacement, h'(s,0,t) to H(s) and express the resulting film thickness, pressure 
and velocities in dimensionless form as 


h = H(S) + h '(S.e.t) , u - 0(S) + u '(S.G.t) , 

v = v(s) + v / (s,e,t) , p = p(S) + p '(s.e.t) . 


(49) 


When above relationships are substituted for the corresponding variables in Equations (40), (41) and (13) and 
quadratic terms in second order (primed) quantities are neglected, the following second order, linearized 
momentum and continuity equations are obtained: 


K / 


1 dp 

fp* 30 


= FT 


/ \ 


du' r.du' dU./ 0 30 7 

+ V + v + 

at 3S dS f 30 


+ tf.u' + tf-v' + <F fi FT , (50) 


r / 


1 3p 1 
p* 3S 


‘ 3v 7 3v ' dV „ / U 3v ' N 
+ V + — v + 


l at 


3S dS f 30 


+ r 0 Q' + T*.v' + , (51) 


vMl + dVfh' + + + + + = o . 

3S dS 3S dS 30 30 at 


S.“ / 




(52) 


The subscripts 0, V and H are used to denote partial differentiation with respect to those variables. 

The continuity requirements for p' and v' at values of s where h is discontinuous are obtained by performing 
a similar linearization on Equation (16) and are given below: 
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v' 


VAh 


H (H - Ah) 


r / ( H - Ah ) _ / 

h + - — : — iv j 

H 


Ap 


*/ - 1 ft 


/. 


^^(XvV^Xnh ) @S = S a 


(53) 


The tangential velocity u' will be continuous. The end conditions given below are obtained by linearizing 
Equations (19) and (20) and substituting the values of U, V and P prescribed by Equations (47) and (48): 

= 0 , p' = -^P*R*(Xv'' / + Xfih ) @s = S in (54) 

and 

P ' = 0 @ S = S ex . (55) 


Solution to first order flow equations 

The first order flow equations are solved by a fairly straight forward numerical procedure. To begin with one 
may integrate Equation (46) and express the result in the form 


fH 

where the dimensionless transverse inlet velocity V im is unknown and must be initially assumed. One may 
compute V(S), from the above equation, for the assumed value of V )n . The tangential velocity, U(S), may now 
be determined by solving Equation (44) numerically subject the first boundary condition in Equation (47). 
Once U(S) and V(S) have been established, P(S) may be determined by numerical integration of Equation 
(45) subject to the second boundary condition in Equation (47). The correct value of V, n will be that which 
results in a pressure distribution that satisfies Equation (48) and is determined numerically by Newton's 
method. The transverse flow rate, Q, may then be calculated as Q = 27ir 0 CV 0 r ln H, n V| n . 

A semi-implicit algorithm has been used for solving Equation (44). Examples of semi-implicit algorithms may 
be found in many texts on numerical methods such as reference 12. Although they are slightly more complex 
than some of the widely used explicit methods such as Runge-Kutta, they provide increased stability for "stiff' 
systems such as that which may occur here at small values of R*. Although Equation (44) is a single first 
order differential equation, the algorithm is presented below for a system of N first order ordinary differential 
equations, as will occur later in the development of the second order solution. 

The system of differential equations may be represented in the form 
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i - 1 ,2,...,N 


(56) 


dY i 

= F i(S.Y v Y 2 Y n ) , 

as 

where Y, denotes the dependent variables. If the values of Y, at the new position (S + AS) are denoted by the 
column vector {Y new }, and the values at the old position by {Y}, the equation used for computing {Y new } may 
be written as 

{Y new> ={Y} + AS([l]-^-[k])-i{F (S+ ^-,Y 1 ,Y 2 Y N ) } , (57) 


where [I] is the unit diagonal matrix of order N and [k] is a matrix of derivatives whose components are 


kij 5Y. 


5 F,(S + ^,Y,,Y 2 Y, Y n ). 


i.j = 1,2 N 


Solution to second order flow equations 


(58) 


Although the global film thickness H, remains an arbitrary function of the transverse position S, the 
displacement of the center of the seal will be taken to have the form of a small oscillatory translation, e x , in 
the x direction, for a cylindrical seal or e z , in the z direction, for a face seal and a small oscillatory rotation 
about the y axis, i|j, for either type of seal. If the seal oscillates at speed Q, the dimensionless film 
displacement, -h', will take the form 


-h ' = |ty o Scos0 + l c e xO cos0 + l f e z0 ]cosQt , 


(59) 


where 


_ Or r 0 

l c = 1 - l f , Q = — -, tJj = — ifr 0 cos^t , e x = — ^-cos&t, e z = — ^ cosOt. (60) 

V 0 C C C 


Each of the three terms on the right hand side of Equation (59) may. be defined as component film 
displacements, 6 k and the equation may be rewritten as 




k=1 


(61) 


The component film displacements may be expressed in complex form as 
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\ = j*(K * *0 . (62) 

where 

. (63) 


The functions 9t(Z) in the above equation and 3(Z) to be used later denote the real and imaginary parts of the 
complex variable Z respectively. The values of e k , q k and J k are given below: 


k = 1. 2, 3 

e k - ^0. ®x0’ ®z0 

r lk = S, l c . I f 

J k = 1, 1, 0 


(64) 


Since the second order flow equations are linear one may superimpose component solutions obtained for h' 
appearing in Equations (50) - (54) replaced by the various values of -6f. If the component solutions for the 
dimensionless second order pressure and velocities are denoted by p£ u k and vj respectively, then the values 
of p', u' and v' may be expressed as 

P / = £Pk- Q/ = £°k- v ; = £v k , (65) 

k=1 k=1 k=1 


where 


Pk = ^(Pk 


P 0 


-9t(u k + + 0 k ) , 


v k = l9t(v k + + v k ) 


( 66 ) 


One may now attempt to separate variables by looking for component solutions in the form 


„ ± 
Pk 


— ± 

Pk 


(S)e' 


(J ‘ e±£S) , 0 k 4 =u k 4 (S)e , ^ ± “ ) 


v* = v k (S)e i(Jk0±ffi) 


(67) 


When the above expressions for p k , u k and v k are substituted for p', u' and v' in Equations (50) - (53) a set 
of ordinary differential equations is obtained which may be arranged in the form 


NAS A/CR— 2003-2 1 2359 


21 



- du k 

R* V — - + 

dS 


<*0 


iR’Q^u* 


- \ 


** + R 


dU 

dS 


— ± 
v k 


i Jl 


P r 


Pk* = e k Tlk^ 


(68) 




1 d(rH)-± 

— ~ v k 

fH CIS 


e k ^ d (~ rT lk) 

rH i cJS 


{ - ) 


+ 

dV . A ± 

"lk 


dS J 


( 69 ) 


1 d P k 
6* ds 



- ) 


r / ~ 


\ 

-I 

¥*.. - iR* — 

— ± 

U k + 

% + R* 

dV 

_ V 

d ( ? H ) +iC f 


l u 

P J 

K 

V 

dS 

fH 

dS 

) 



v k 


= e * r n\ ~ e k R 


♦ _v_ 

rH 


+ 

dS 


dV . A ± 
— + iQ 

^ dS 


H k 


where 


Q* = ±Q + — 


( 70 ) 


( 71 ) 


The transformations given in Equation (66) satisfy the periodicity requirements in 6. The remaining boundary 
and continuity conditions are obtained by substituting the variables defined in Equation (67) for the 
corresponding variables in Equations (53) - (55). The form of the equations remains unchanged and the 
results are given below: 


v t = — 


VAh 


H (H - Ah) 


— ' Vlk + 


(H ". Ah> v t * J . Ap„* = )p*R , (XvV l *-x fi e,,%) @s = S, , ( 72 ) 

H 2 


U k = 0 . Pk = ~P*R*(Xv V k - Xh 6 ^) @ S = S in 
and 

Pk = o @ S = s ex . 


( 73 ) 

( 74 ) 


Equations (68) - (70) are solved numerically using the algorithm defined by Equations (56) - (58). 
Complimentary and particular solutions are constructed and combined so that the boundary conditions given 
by Equations (72) - (74) can be satisfied without having to iterate. 
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The real component pressures p k , can be expressed in terms of the real and imaginary parts of p k and p k by 
expanding the complex pressures, p£, substituting the result for the complex pressures in Equation (66) and 
extracting the real part as indicated. The result may be written as 

Pk = "[^(Pk* + Pk )cosJ k 0 - S(p k * + p k )sinJ k 0]cosQt - j[s(p k + - p k ")cosJ k 0 + 3t(p k * - p k ~)sinJ k 0]sinQt . (75) 

The first term on the right in the above equation will be in phase with the displacement and related to the 
stiffness of the seal. The negative of the second term will be in phase with oscillatory velocity and related to 
the damping coefficient. 

Forces, moments and dynamic coefficients 

The dimensionless applied loads and moments obtained by integrating the pressure over the seal area are 

2 It S R 2jt S R 2 It S R 

w x = l c J J pcos0rdSd0, w y = l c J J psin0rdSd0, w z = l f J J prdSd0, 

0 s L 0 s L 0 S L 

(76) 

2ti §r 2n s r 

M x = - j" J psin0rSdSd0, M y = J J pcos0rSdSd0 . 

0 S L 0 S L 

All first order applied loads and moments will be zero with the exception of the first order axial load acting on 
a face seal, W z , which is given in dimensionless form by 

s R 

W z = 2 tcI f f PrdS . (77) 

s L 

The various stiffnesses can be obtained by substituting the first term on the right hand side of Equation (75) 
for the pressure, p, appearing in the force and moment integrals in Equation (76) and dividing the result by 
the associated displacement of the center of the seal. The dimensionless translations and rotations that have 
been included in Equation (59) are a rotation about the y axis of magnitude ij; 0 cosQt (corresponding to k=1), 
a translation in the x direction of magnitude e x0 cosQt (corresponding to k=2) for a cylindrical seal and a 
translation in the axial direction of magnitude e z0 cosQt (corresponding to k=3) for a face seal. The damping 
coefficients may be obtained in a similar manner by substituting the second term on the right hand side of 
Equation (75) for the pressure in Equation (76) and dividing the result by the velocity -Qt ijj 0 sinQt for k=1, 
-Qe x0 sinQt for k=2 and -Qe z0 sinQt for k=3. 

For small displacements about the centered position, stiffness and damping values associated with 
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translations in the y direction or rotations about the x axis may be obtained from symmetry considerations. 
For example, a displacement in the y direction will produce a force in that direction of the same magnitude 
as that which would be produced in the x direction for an equivalent displacement in that direction (K yy =K xx ). 
Similarly, a displacement in the y direction will produce a force at 90° from that direction (which would be the 
negative x direction) as would be produced in the positive y direction for an equivalent displacement in the 
x direction (K xy =-K yx ). 

The symbols <j) and i|r, when appearing as second subscripts on the dynamic coefficients, will denote rotations 
about the x and y axes respectively. When appearing as first subscripts they will denote moments about the 
x and y axes. The dimensionless stiffness and damping coefficients relating moments about the x and y axes 
to a rotation about the y axis (k=1) for both face and cylindrical seals are 


K It 

= »(P,* + P,')SdS , K „ = -£-/ 8t<p,* + p,-) 

C- 1 4 


SdS , 


s„ 


(78) 


B 


J $ ( p,*-p, >sds . B w = -^-|3(p,*-pOS 


dS . 


1 S, 


s L 


The corresponding coefficients for a rotation about the x axis obtained from symmetry considerations are 




(79) 


The dimensionless dynamic coefficients relating the x and y force components to a rotation about the y axis 
(k=1) for a cylindrical seal are 


^xv|; 


K 

2e„ 


f + pO 


dS , K y(|f = - 


s L 


2^ 


f ^(p; + Pl )ds 


B 


Sr S 

7Cl . _ - ItL 

c r rs>,_ + -. \ jo o _ c 


Xl|; 




f S(p,*-p,-)dS . B rt 


2Qe 


f 9t( P ;- Pl -)ds 


i s. 


(80) 


and the corresponding coefficients for a rotation about the x axis are 
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K xO> ^y»|j 1 ^y<j> ^xv|i ' ®x(J) ®yi|; ' ®y4> 


-B 


xi|; 


( 81 ) 


The dimensionless stiffness and damping coefficients relating moments about the x and y axes to a translation 
in the x direction (k=2) for a cylindrical seal are 


pi pi 

K = 3(Pz *P 2 ‘)SdS , = / »(P 2 * + P 2 ') 

62 S L 62 S L 


SdS , 


B 


4>x 


pi 

7— f ®(P2 -P 2 ") 
^ J 


2&e 


SdS ■ = f S(P 2 *-P 2 ‘) 

2Qe J 


SdS 


2 S L 


2 S. 


(82) 


and the corresponding coefficients for a translation in the y direction are 


^^y K.X , K^y , B^ B^ , B^y B^ 


(83) 


The dimensionless dynamic coefficients relating the x and y force components to a translation in the x 
direction (k=2) for a cylindrical seal are 


Kxx = 


B„„ = 


2e„ 


n\. 


iz\. 


f 9t( P ; + p 2 )dS , Ky, = — £- I 8(p 2 + + p 2 )dS 


Til 


2&6 


f S( p 2 + -p 2 )ds, § = — S- f gt(p 2 + -p 2 )ds 

J 2Qe J 


2 S. 


2 S, 


(84) 


and the corresponding coefficients for a translation in the y direction are 


K = — K k — K R — _d p — p 

^xy ^yx ' ^yy ^xx 1 D xy D yx ’ D yy u xx 


(85) 


Finally, the stiffness and damping coefficients relating the axial force to an. axial translation (k=3) for a face 
seal are 
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( 86 ) 


TC I R j L | ° 

R zz = — f Sft(p 3 + +P 3 ')dS , B zz = —t- f ^(p 3 + -p 3 “)dS . 

63 S L ^ 6 3 S L 

As a result of the assumption of small displacements about the centered position used in this analysis, the 
remaining dynamic coefficients for a face seal will all be zero: 

K,* = K*. = K* = K* = B z4 , = = B z „ * B t; = 0 . (87) 

Relationships of the following type may be used to calculate the physical stiffness and damping coefficients 
from the dimensionless ones 


and 


K _ P 0 r 0 £ v 

^XX ^ Nx ' q 


Po r o " 


K, 


K 


X<j> 


3 

>o r o 


K 


xc|> 


B„ 


Po r o 

CV„ 


B.. 


Bj 


Po r o 

CV„ 


B 


<t>4> 


B 


XtJ> 


Po r o 

CV„ 


B 


X<J) 


(88) 


(89) 


The stiffness coefficients, which relate to forces in phase with the displacements and the damping coefficients, 
which relate to forces 90° out of phase with the displacements, are both, in general, functions of Q. For the 
case of incompressible flow under consideration, the damping coefficients and the cross coupled stiffness 
coefficients (K^, K^, Ky X , K^) are nearly independent of Q. Under these circumstances, it has been found 
that the frequency dependence may be described adequately with the introduction of an apparent mass matrix 
[A] such that [K] = [K°l - Q 2 [A], where [K°] is the stiffness matrix evaluated at £1=0. The components of [A] 
may be evaluated by computing [K] at a characteristic speed such as Q=V„/r 0 corresponding to Q=1 or 
synchronous speed (Q=to) and using the formula 


K - K° 

A„ * -a '-L (90) 

J Q 2 

The cross coupled components of [A] which have been included for completeness are not significant and the 
direct components are indeed found to be nearly independent of £1. 

Although the linearizations that have been used in obtaining the second order flow equations are only valid 
for small values of e k , once they have been performed the resulting equations for the dynamic coefficients may 
be shown, by rescaling variables, to be independent of the value e k actually used. Any convenient value such 
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as e k =1 may be used in computations. 


Determination of power loss and flow 

When the rotating surface (a) is ungrooved, the power loss will be proportional to the integral of rx a -T over 
the surface of the seal. When the rotating surface is grooved, an additional term will be needed to correct for 
the pressures acting edges of the grooves. Although the solution to the equations based on bulk flow theory 
used here does provide for the prediction the shear stress on each surface, the validity of power loss 
predictions is somewhat doubtful. This is due to the difference in characterizing the bulk flow behavior for 
velocity distributions associated with Couette and Poiseuille flow. 

One may write x a as % = j + p x wherq. x denotes the Couette part of the shear stress defined as 
x c = (x a + x b )/2 and x p denotes the Poiseuille portion of the shear stress defined as x p = (x a - x b )/2. It can 
be seen from Equations (1) and (2) that only x p appears in the primary and secondary flow equations used 
to predict load, flow and dynamic coefficients. Comparisons between theory and experiments based on 
measurements of flow and dynamic coefficients alone may be used to validate or determine values of x p , but 
will not necessarily result in a bulk flow model that accurately predicts power loss. This weakness in the bulk 
flow model was recognized by Hirs (reference 1) who suggested the use of different friction factor 
representations based on measurements obtained from experiments based on Couette and Poiseuille flow. 
The method for predicting power loss used here will be a highly approximate one based on generalizations 
of an approach used by Hirs for cases where both surfaces were smooth and the flow was predominantly in 
the tangential direction. 

If the scaler shear stress, x, is used to represent the tangential shear stress, (x=x-T), then the dimensionless 
forms of the tangential shear stresses x a and x b , can be obtained from Equations (4), (5), (8), and (10) as 

% = , % = -P'RA(R b )T • (91) 

h h 


One may attempt to compensate for the influence of the Couette portion of the shear stress by replacing x a 
with the effective shear stress, x*, defined as 


t_ + 


2 A 



(92) 


in the integral used for determining power loss. The correct value of the "Couette reduction factor", k, may 
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be shown to be A.=3 for laminar flow. A value of A,=1.2 will be used for turbulent flow. This value was 
suggested by Hirs based on one dimensional Couette flow measurements and provides power predictions 
that are in good agreement with those obtained with Ng and Pan theory (reference 13) over its range of 
validity. 

The dimensionless power loss, P = P/(Cp 0 r§cD), obtained from the solution to the first order flow equations is 
given by 



The quantities x g and Xr appearing in the above equation denote the values of the effective shear stress, x*. 
based on the first order velocities and film thicknesses for the grooves and ridges respectively. The effective 
shear stress is defined by Equations (91) and (92) and the local velocities are obtained from the solution to 
Equations (30) and (32) - (34) and the definitions provided by Equation (21). The sum of the first two terms 
appearing inside the brackets of the integrand in Equation (93) represent the global value of x*. The third term 
has been added to include the forces acting on the edges of the grooves when the grooves are on the rotor. 
The quantity P g e is obtained by substituting first order variables for the corresponding quantities in Equation 
(24) and is given below: 


9.0 


= -P** g - P*R* 


v^ + ,. uv 


dS 


■r - 

r 


( 94 ) 
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3. Description of Computer Code SPIRALI 


The FORTRAN computer program SPIRALI, has been written to implement the analysis that is contained in 
Section 2. The program has been compiled in its present form with Version 5.1 of the Microsoft® Fortran 
Compiler and should work with many other compilers with relatively little modification. SPIRALI uses two 
features that are supported by a wide variety of compilers that are extensions to FORTRAN 77. The first of 
these extensions is the use of NAMELIST input. The second extension relates to the implementation of 
complex double precision variables in the solution of the second order flow equations required to compute 
stiffness, damping and apparent mass values. It should be noted that NAMELIST is a VAX extension and 
COMPLEX DOUBLE PRECISION is part of IBM Systems Application Architecture, IBM VS and VAX 
extensions to Fortran 77. 

The main program is system dependent in that it opens files and when compiled with the Microsoft compiler, 
allows placement of file names on the command line. It also reads the input data and generates the plot file. 
The supporting output routines either write to the standard output unit or take FORTRAN supported logical 
unit numbers as input, hence should not require significant modification with any FORTRAN 77 compiler. In 
particular, the subroutine INLIST writes the NAMELIST to the output file in a legible manner and the remainder 
of the output file is generated by the subroutine DIMOUT. Status and error messages are written to the 
standard output unit by the subroutines OUTSCR and EMSG respectively. There are no input or output 
statements in any of the other routines contained in SPIRALI. 

Although SPIRALI is designed to run over a full range of conditions of laminar or turbulent flow, certain 
precautions need to be taken when the transverse flow rate is small. It can be seen, for example, that 
Equation (44) becomes singular as R'V approaches 0. Numerical problems can be avoided here by setting 
RV=0 and numerically solving Equation (44) as an algebraic equation rather than as a differential equation. 
SPIRALI provides an option for the user to exclude the transverse inertia effects arising from the terms v3u/3S 
and l f uv/r in Equation (11) and v3v/3S in Equation (12) and the associated discontinuities in pressure that 
occur at the inlet and at jump discontinuities in film thickness. This option, which is controlled by the input 
parameter NOI (see input description) is particularly useful for situations where the rotation speed is high but 
the seal pressure difference is low. 

If transverse inertia is included, the location of the upstream or inlet boundary must be specified. When the 
seal pressure difference is sufficiently high, the upstream boundary will be at the high pressure side of the 
seal. At lower pressure differences the upstream boundary may depend on the rotor speed and the spiral 
groove parameters. The input parameter IFLOW is used to select the upstream boundary. If you are not sure 
of the flow direction, run the program with transverse inertia excluded. If you then wish to include transverse 
inertia, adjust IFLOW in accordance with the sign of the transverse flow rate given in the output and rerun the 
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program with transverse inertia included. 


The analytical procedure contained in Section 2 has been oriented toward determining pressure distribution, 
load, flow, power loss, stiffness and damping for a given centered film thickness distribution. Face seals will, 
in general, have a non-zero first order axial load component and it is often desirable to determine the 
equilibrium film thickness from a prescribed load used to balance the seal. SPIRALI provides an iterative 
homing option for this purpose that is controlled by the input parameter IHOME. The homing option should 
be used with caution since convergence is sensitive to the starting film thickness which is determined from 
the input parameter C. It is recommended that the program be run for various values of C without the homing 
option so that a good starting value of C (one that produces an output value of the load to balance the face 
seal that is reasonably close to the prescribed load) can be obtained for the homing process. 



Figure 7 Parameters for characterizing quadratic film variation. 


Although the analysis provides for an arbitrary transverse film thickness variation, SPIRALI is limited to a 
quadratic film variation with coefficients determined from the parameters H TAP , H BRL and C as shown in Figure 
7. These parameters may be positive, as shown or negative and should be able to characterize a fairly 
general film variation without the complexity of required inputs at every grid point. 
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Print inputs and check 
for illegal values. Any 
errors? 


Generate grid. 

I 

Compute film and trans- 
verse velocity scale fac- 
tor. Any errors? 

t n 

Compute dimensionless 
variables and parame- 
ters. Any errors? 

- - |l; - : i 

Set flow direction. 


Solve first order flow 
equations from inlet to 
exit. Any errors? 


Set new inlet velocity. 


Is exit pressure con- 
dition satisfied? 

u ~ 

Compute load, flow 
and power loss. 


Has iteration limit 
been reached? 


n Compute new value 
— of C. Are C and axial 
stiffness positive? 


Has iteration limit 
been reached? 



Solve second order flow 
equations at 0 frequency 
and at running freqency. 
Any errors? 


Compute stiffness, 
damping and mass 
matrices. 


Home in on axial load for 
face seals? 


Print available output 
and exit. 


Figure 8 Flow diagram for overall logic used in computations 
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The logic used in SPIRALI for solving the first and second order flow equations, computing load, flow, power 
loss, stiffness, damping, apparent mass coefficients, etc. for a given case is shown in Figure 8. On 
completion, the program will look for a new NAMELIST and if found, will proceed to the next case. 

The program SPIRALI may be run in its present form under either DOS or OS2 operating systems by entering 
the command 

SPIRALI infile outfile plotfile 

where infile and outfile are user supplied names of the input and output files respectively and plotfile is a data 
file containing first order film thickness, velocity and pressure profiles as functions of S. The contents of these 
files are described below. 


Input description for SPIRALI 

SPIRALI requires its input in NAMELIST form. The NAMELIST variables will be given below in the order that 
they are listed in the output. Units for dimensional inputs will be given in English units as required when 
ISIUN=0, followed by SI units required when ISIUN=1 in parentheses, where applicable. The NAMELIST 
name is INPUTS. The relationship between the NAMELIST variables and symbols appearing in the 
nomenclature will be given wherever possible. 


TITLE 
I FACE 
ISIUN 
I GROT 
NOI 


I FLOW 


RO 

EL 


any character string enclosed in quotes containing up to 64 characters 
0 for cylindrical seal, 1 for face seal 
0 for English units, 1 for SI units 

groove rotation flag, set to 1 , if the rotating surface is grooved or 0, if ungrooved 
parameter which controls treatment of inertia effects with following values: 

-1 include all inertia effects 

0 include tangential inertia effects, program decides on transverse inertia effects 

1 include tangential inertia effects, exclude transverse inertia effects 

2 exclude all inertia effects 

parameter which controls selection of inlet boundary with following values: 

-1 assume transverse flow rate to be negative, inlet boundary is at right (s in =s R ) 

0 assume that the sign of the transverse flow rate is the same as the sign of p L -p R and 
set inlet boundary accordingly 

1 assume transverse flow rate to be positive, inlet boundary is at left (S|„=s L ) 
r 0 , radius of cylindrical seal (IFACE=0) or outside radius of face seal (IFACE=1), in. (m) 

L, length of cylindrical seal (IFACE=0) or land width of face seal (IFACE=1) as shown in Figure 
1, in. (m) 
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c 


RPM 

RPMO 

RPMD 


PLEG 

PRIG 

FZD 


vise 

DENS 

EMA 

ENA 

EMB 

ENB 

HTAP 

HBRL 

TOLH 

TOLV 

DUT 

IHOME 


NITH 

NITV 

NREG 

NRSUB 

ELFR 


C, constant spacing between the two surfaces of the seal which would occur if 6, H TAP and H BRL 

were all equal to 0 (see Fig. 7). and is referred to as the clearance of a cylindrical seal (IFACE=0) 

or nominal film thickness of a face seal (IFACE=1), in. (m) 

w, angular speed of rotating surface converted to revolutions per minute 

u,„/r ln , angular swirl speed of fluid at inlet to seal converted to revolutions per minute 

Q, angular speed of disturbance converted to revolutions per minute. If Q>0, the output file will 

contain a table of stiffness and damping values, evaluated at Q followed by the static stiffnesses 

evaluated at 0 frequency. If Q=0, the static stiffnesses will be replaced by apparent masses. 

p L , gage pressure at s=s L , psi (Pa) 

p R , gage pressure at s=s R , psi (Pa) 

W z , applied load in z direction, used only if IH0ME=1 and IFACE=1, lb. (N). It should be noted 
that the ambient pressure, min(p L ,p R ), does not contribute to the applied load and should be 
subtracted from all pressures when determining pressure induced loading, 
p, fluid viscosity, psi-sec (Pa-sec) 
p, fluid density, lb-sec 2 /in 4 (Kg/m 3 ) 

m a , exponent in shear factor power law relationship for moving surface 

m a , coefficient in shear factor power law relationship for moving surface 

m b , exponent in shear factor power law relationship for stationary surface 
m b , coefficient in shear factor power law relationship for stationary surface 
H tap , film change defining linear taper, see Fig. 7, in. (m) 

H brl , film change defining quadratic "barreling", see Fig. 7, in. (m) 

maximum allowable relative error in axial load (used only when IFACE=1 and IHOME=1) an error 
code is returned if this tolerance is not met in NITH iterations 

maximum allowable relative error in dimensionless first order transverse inlet velocity, V, n . An 
error code is returned if this tolerance is not met in NITV iterations, 
relative increment used in numerical differentiation 

when set to 1 triggers computation of the nominal film thickness C for which the applied load 
prescribed by FZD is equal to the calculated load required to balance the face seal (used only 
if IFACE=1) 

maximum number of iterations to be used for homing in on load for face seal (used only if 
IFACE=1 and IHOME=1) 

maximum number of iterations to be used determination of first order transverse velocity 
component at inlet, V in 

number of regions in axial direction (IFACE=0) or radial direction (IFACE=1) 
vector of length NREG containing number of sub-intervals in each of the NREG regions 
vector of length NREG containing fractional extent of each of the NREG regions (length of region 
divided by L). The sum of the NREG values of ELFR must add up to 1 . 


NAS A/CR— 2003-2 1 2359 


33 



ALPI vector of length NREG containing groove to pitch ratio, a, for each of the NREG regions. Set a=0 
in ungrooved regions. Set a=1 for regions comprising circular steps or grooves. Spiral groove 
theory is not implemented when a=0 or a=1. 

BETI vector of length NREG containing groove angle, (i (degrees), in the range -90° < BETI < 90°, for 
each of the NREG regions. Set to 0 in ungrooved regions. Spiral groove theory is not 
implemented when p=0. 

DELT vector of length NREG containing groove depth, 6, for each of the NREG regions, in (m). If a=1 
but 6*0, the region will be treated as a circular groove (6>0) or step (6<0). 

ZET vector of length NREG containing the contraction loss coefficient, C, at the inlet to each of the 
NREG regions 

NSG vector of length NREG containing the number of spiral grooves in each region. This input is used 
solely for characterizing effects of local inertia on pressure jumps in regions where spiral groove 
theory is implemented. In regions where NG=0 effects of local inertia on spiral groove 
performance is ignored. 

ZETG vector of length NREG containing the contraction loss coefficient, £, associated with the local 
pressure drop for flow going from grooves to ridges for each of the regions where spiral groove 
theory is implemented and NG>0. 

The default values of the NAMELIST variables are given in Figure 9. These values correspond to a static, 
ungrooved, cylindrical seal with inputs in English units. The default case cannot be run until non-zero values 
are provided, at least, for R0, EL, C, VISC and either PLEG, PRIG or RPM. 

Description of output from SPIRALI 

Listings of output files generated by SPIRALI are given in Section 4 under sample problems. The output file 
starts with a reference case number and the title which is followed by a listing of the NAMELIST variables in 
a form suitable for extraction as an input file. The next line gives the seal type and a statement as to which 
(if any) inertia effects have been excluded. The next 4 lines provide a brief description of the seal geometry, 
operating data and fluid properties. When IHOME=0 and IFACE=1 the nominal film thickness will be the 
computed value rather than the input value of C. The next line contains an error code (see Section 7 for an 
explanation of non-zero error codes) and the number of iterations used in determining the transverse flow rate 
when solving the first order flow equations (iterations in primary flow). When IFACE=1 the axial load required 
to balance the face seal will be given which may be compared with the input value of FZD, when IHOME=1 , 
to evaluate the accuracy of the homing process. The next two lines contain the transverse flow rate, Q, the 
torque and the power loss P. These are followed by the local Reynolds numbers at the left and right seal 
boundaries based on the absolute values of the transverse fluid velocities and the circumferential fluid 
velocities relative to the rotor respectively and may be used to estimate the degree of turbulence in the seal. 
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&INPUTS 


TITLE 

= 

i i 






I FACE 

= 

0 

ISIUN 

= 

0 



I GROT 

= 

0 

NOI 

= 

0 

I FLOW 

= 0 

RO 

= 

0. 0000E+00 

EL 

= 

O.OOOOE+OO 

C 

= 0.0000E+00 

RPM 

= 

0.0000E+00 

RPMO 

= 

O.OOOOE+OO 

RPMD 

= 0.0000E+00 

PLEG 

= 

0.0000E+00 

PRIG 

= 

O.OOOOE+OO 

FZD 

= 0.0000E+00 

vise 

= 

0. 0000E+00 

DENS 

= 

O.OOOOE+OO 



EMA 

= 

-2.5000E-01 

ENA 

= 

7.9100E-02 



EMB 

= 

-2.5000E-01 

ENB 

= 

7.9100E-02 



HTAP 

= 

O.OOOOE+OO 

HBRL 

= 

O.OOOOE+OO 



TOLH 

= 

1.0000E-04 

TOLV 

= 

1 . 0000E-05 

DUT 

= 1.0000E-06 

I HOME 

= 

0 

NITH 

= 

10 

NITV 

= 30 

NREG 

= 

1 

NRSUB 

= 

20 



ELFR 

= 

1.0000E+00 






ZET 

= 

O.OOOOE+OO 






ALPI 

= 

0.0000E+00 






BETI 

= 

O.OOOOE+OO 






DELT 

= 

O.OOOOE+OO 






NSG 

= 

0 






ZETG 

= 

O.OOOOE+OO 







/ 


Figure 9 Namelist defaults 


A table of stiffness, K, damping, B, and apparent mass A, (if RPMD=0) or static stiffness, KO (if RPMD>0) are 
given next, where the column corresponds to the displacement or rotation and the row corresponds to the 
force or moment. The units for any given stiffness, damping or mass value will be the force unit given at the 
end of its row divided by the displacement unit given at the top of its column. 


The plot file corresponding to the third argument on the command line starts with a line containing the number 
of points in the data set to be plotted, NP. This is followed by a table of NP lines. The columns from left to 
right are values of the variables S, H r [in. (m)], U [in/sec (m/sec)], V [in/sec (m/sec)] and P [psi (Pa)]. 
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4. Sample Problems 


Several sample problems have been prepared to show the use of the computer program for cylindrical and 
face seal applications. They are intended primarily for illustration and do not necessarily represent 
recommended seal designs. Separate input files are not given since the complete input NAMELIST is 
included in the output file for each case. 

Cases 1 - 5 apply to the cylindrical seal geometry shown schematically in Figure 10. The seals for cases 1 - 
4 are divided into two regions of equal length. The region at the left consists contains a spiral groove pattern 
for cases 1 and 2 and an ungrooved step of depth 6 for cases 3 and 4. A plain cylindrical seal consisting of 
a single region at a constant film thickness, C, is used for comparison in case 5. The high pressure is at the 
left hand side of the seal for cases 1 , 3 and 5 and at the right hand side for cases 2 and 4. A groove angle 
of p=25° has been used for cases 1 and 2, which will tend to pump from left to right. The flow produced by 
the pressure gradient should be in the same direction as groove pumping in case 1 and in the opposite 
direction in case 2. The pressure gradient will tend to dominate the flow in both cases. The parameters NOI 
and IFLOW have both been set to 0 for all of the above cases. Inertia effects have been included for all of 
the sample cases (otherwise the output file would state which effects were neglected) and the assumption that 
the flow direction is away from the high pressure side of the seal, which is employed when IFLOW=0, is 
validated by the sign of the flow rate for each of the cases. 

Cases 6-8 correspond to a mechanical face seal with spiral grooves on the outside surface of the stator 
oriented to pump inward. The stator geometry is shown schematically in Figure 1 1 . The high pressure is on 
the outside of the seal. A 1 mil (.0254 mm) nominal film thickness is used for case 6 and a value of 1.5 mils 
(.0381 mm) is used for case 7. IHOME is set equal to 1 for case 8 and a load of 1600 lb (71 16.8 N) is applied 
to balance the seal. The load to balance the seal used in case 8 was bounded by the results of cases 6 and 
7. This bounding procedure is recommended for all applications where IHOME will be set to 1 . A finer grid 
was used in the inlet (outside) region for cases 6 - 8, to reduce truncation error resulting from the inclusion 
of inertia effects with the relatively low radial Reynolds number that occurs at the outside radius in case 6. 
It is particularly important to check accuracy by increasing the values of NRSUB when inertia effects are 
included. 

Cases 9 - 10 correspond to a straight through helically grooved cylindrical seal of the type reported by Childs, 
Nolan and Kilgore in Reference 14. The input geometry for Case 9 corresponds to their Stator #2 which has 
17 grooves, a groove angle of 15°and a gap of .376 mm. Case 10 is the same as Case 9 except that the 
value of N g has been set equal to 0. Case 10 does not have any real physical significance and is presented 
solely for comparison with Case 9. 
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Figure 10 Schematic of cylindrical seal for cases 1-5 
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( CASE 1 ) Cylindrical seal, grooves pumping inward, high pres, at left 
&INPUTS 

TITLE = 'Cylindrical seal, grooves pumping inward, high pres, at left' 


I FACE 

= 

0 

ISIUN = 

0 



I GROT 

= 

0 

NOI 

0 

I FLOW = 

0 

R0 

= 

1 . 0000E+00 

EL 

5 . 0000E-01 

C 

1 . 000 0E- 03 

RPM 

= 

5 . 0000E+04 

RPM0 

2 . 5000E+04 

RPMD 

0 . 0000E+00 

PLEG 

= 

1 . 0000E+03 

PRIG 

0 . 0000E+00 

FZD 

0 . 0000E+00 

vise 

= 

3 . 0000E-08 

DENS 

1. 0000E-04 



EMA 

= 

- 2 . 5000E- 01 

ENA 

7. 9100E-02 



EMB 

= 

- 2 . 5000E- 01 

ENB 

7 . 9100E- 02 



HTAP 

= 

0 . 0000E+00 

HBRL 

0 . 0000E+00 



TOLH 

= 

1.0000E-04 

TOLV 

1 . 0000E-05 

DUT 

1.0000E-06 

I HOME 

= 

0 

NITH 

10 

NITV 

30 

NREG 

= 

2 

NRSUB = 

50 50 



ELFR 

= 

5 . 0000E-01 

5.0000E-01 




ZET 

= 

0 . 0000E+00 

0 . 0000E+00 




ALP I 

= 

5.0000E-01 

0 . 0000E+00 




BETI 

= 

2 . 5000E+01 

0 . 0000E+00 




DELT 

= 

2 . 0000E-03 

0 . 0000E+00 




NSG 

= 

16 0 





ZETG 

= 

0 . 0000E+00 

0 . 0000E+00 





/ 

CYLINDRICAL SEAL 

LENGTH, DIAMETER, CLEARANCE = 5.0000E-01, 2.0000E+00, 1.0000E-03 (IN) 

ROTOR, SWIRL AND DIST. SPEEDS = 5.0000E+04, 2.5000E+04, 0.0000E+00 (RPM) 

PRESSURE AT START, END AXIAL BOUNDARIES = 1.0000E+03, 0.0000E+00 (PSI) 

VISCOSITY = 3.0000E-08 (PSI-SEC) , DENSITY = 1.0000E-04 (LB-SEC/IN4) 

ERROR CODE = 0, ITERATIONS IN PRIMARY FLOW = 5 

FLOW = 9.6718E+00 (IN**3/SEC) 

TORQUE = 6 . 2612E+00 (IN-LB), FILM POWER LOSS = 4.9672E+00 (HP) 

AXIAL REYNOLDS NUMBER = 1.0262E+04 

CIRC. REYNOLDS NUMBERS FOR ROTOR AT SEAL ENDS = 3.4907E+04, 1.7454E+04 

DYNAMIC COEFFICIENTS ( FORCE UNIT / DISP. UNIT ) 


DISP. 

x (IN) 

y (IN) 

phi (RAD) 

psi (RAD) 

FORCE UNIT 

Kx 

3 . 6100E+05 

1 . 9342E+05 

-1 . 6886E+04 

1 . 2744E+05 

LB 

Ky 

-1 . 9342E+05 

3 . 6100E+05 

- 1 . 2744E+05 

-1 . 6886E+04 

LB 

Kphi 

-9 . 9981E+02 

2 . 1794E+03 

3 . 1907E+03 

8 . 8432E+02 

IN- LB 

Kpsi 

-2 . 1794E+03 

- 9 . 9981E+02 

-8 . 8432E+02 

3 . 1907E+03 

IN- LB 

Bx 

6 . 8435E+01 

1 . 2504E+01 

8 . 1569E- 03 

5 . 6432E+00 

LB- SEC 

By 

- 1 . 2504E+01 

6 . 8435E+01 

- 5 . 6432E+00 

8.1569E-03 

LB-SEC 

Bphi 

3 . 8905E- 02 

- 1 . 4517E- 01 

4 . 4703E- 01 

6 . 3076E-02 

IN-LB-SEC 

Bpsi 

1 .4517E-01 

3 . 8905E- 02 

-6 . 3076E- 02 

4 . 4 703E- 01 

IN- LB -SEC 

Ax 

2 . 2877E- 03 

- 7 . 4684E- 05 

- 1 . 3827E- 05 

8 . 266 8E- 06 

LB-SEC2 

Ay 

7 . 4684E- 05 

2 . 2877E- 03 

- 8 . 2668E- 06 

- 1 . 3 827E- 05 

LB-SEC2 

Aphi 

8 . 8069E- 07 

-2 . 7002E- 05 

9. 5112E-06 

-4 . 5208E- 07 

IN-LB-SEC2 

Apsi 

2 . 7002E- 05 

8 . 8069E- 07 

4.5208E-07 

9 . 5112E- 06 

IN-LB-SEC2 


NASA/CR— 2003-212359 


38 



( CASE 2 ) Cylindrical seal, grooves pumping inward, high pres, at right 
&INPUTS 

TITLE = 'Cylindrical seal, grooves pumping inward, high pres, at right' 


I FACE 

= 

0 

ISIUN = 

0 



I GROT 

= 

0 

NO I 

0 

I FLOW 

= 0 

R0 

= 

1. 0000E+00 

EL 

5 . 0000E- 01 

C 

= 1.0000E-03 

RPM 

= 

5 . 0000E+04 

RPM0 

2 . 5000E+04 

RPMD 

= 0 . 0000E+00 

PLEG 

= 

0 . 0000E+00 

PRIG 

1 . 0000E+03 

FZD 

= 0 . 0000E+00 

vise 

= 

3.0000E-08 

DENS 

1. 0000E-04 



EMA 

= 

-2.5000E-01 

ENA 

7 . 9100E- 02 



EMB 

= 

-2 . 5000E- 01 

ENB 

7 . 9100E- 02 



HTAP 

= 

0.0000E+00 

HBRL 

0 . 0000E+00 



TOLH 

= 

1.0000E-04 

TOLV 

1. 0000E-05 

DUT 

= 1.0000E-06 

I HOME 

= 

0 

NITH 

10 

NITV 

= 30 

NREG 

= 

2 

NRSUB = 

50 50 



ELFR 

= 

5 . 0000E-01 

5.0000E-01 




ZET 

= 

0 . 0000E+00 

0 . 0000E+00 




ALPI 

= 

5 . 0000E- 01 

0 . 0000E+00 




BETI 

= 

2 . 5000E+01 

0 . 0000E+00 




DELT 

= 

2.0000E-03 

0 . 0000E+00 




NSG 

= 

16 0 





ZETG 

= 

0 . 0000E+00 

0 . 0000E+00 





/ 

CYLINDRICAL SEAL 

LENGTH, DIAMETER, CLEARANCE = 5.0000E-01, 2.0000E+00, 1.0000E-03 (IN) 

ROTOR, SWIRL AND DIST. SPEEDS = 5.0000E+04, 2.5000E+04, 0.0000E+00 (RPM) 

PRESSURE AT START, END AXIAL BOUNDARIES = 0.0000E+00, 1.0000E+03 (PSI) 

VISCOSITY = 3.0000E-08 (PSI-SEC) , DENSITY = 1.0000E-04 (LB-SEC/IN4) 

ERROR CODE = 0, ITERATIONS IN PRIMARY FLOW = 5 

FLOW = -5 . 8701E+00 (IN**3/SEC) 

TORQUE = 7 . 0124E+00 (IN-LB), FILM POWER LOSS = 5.5631E+00 (HP) 

AXIAL REYNOLDS NUMBER = 6.2283E+03 

CIRC. REYNOLDS NUMBERS FOR ROTOR AT SEAL ENDS = 4.6032E+04, 1.7453E+04 

DYNAMIC COEFFICIENTS ( FORCE UNIT / DISP. UNIT ) 


DISP. 

x (IN) 

y (in) 

phi (RAD) 

psi (RAD) 

FORCE UNIT 

Kx 

1 . 8559E+05 

1 . 6560E+05 

9 . 2632E+02 

-1 . 2796E+05 

LB 

Ky 

-1 . 6560E+05 

1 . 8559E+05 

1 . 2796E+05 

9 . 2632E+02 

LB 

Kphi 

4 . 0924E+03 

- 9 . 6714E+03 

-1.7721E+03 

1. 0930E+03 

IN- LB 

Kpsi 

9 . 6714E+03 

4 . 0924E+03 

-1. 0930E+03 

-1.7721E+03 

IN- LB 

Bx 

5. 9485E+01 

1 . 03 99E+01 

-2 . 8343E- 01 

-2 . 8982E+00 

LB -SEC 

By 

-1. 03 99E+01 

5 . 9485E+01 

2 . 8982E+00 

-2 . 8343E-01 

LB-SEC 

Bphi 

-4 . 013 0E- 02 

-1. 5876E+00 

3 . 1133E- 01 

3 . 9026E- 02 

IN-LB-SEC 

Bpsi 

1. 5876E+00 

-4 . 013 0E- 02 

-3 . 9026E- 02 

3 . 1133E- 01 

IN-LB-SEC 

Ax 

2 . 5027E- 03 

- 8 . 4394E- 05 

6 . 3797E-06 

9 . 3 042E- 05 

LB-SEC2 

Ay 

8 . 43 94E- 05 

2 . 5027E-03 

-9.3 042E- 05 

6 . 3797E- 06 

LB-SEC2 

Aphi 

5 . 0318E-06 

- 1 . 6707E- 05 

9 . 1049E- 06 

3 . 0227E- 07 

IN-LB-SEC2 

Apsi 

1.67 0 7E- 05 

5 . 0318E- 06 

-3 . 0227E- 07 

9.1049E-06 

IN-LB-SEC2 
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( CASE 3 

) 

Cylindrical 

step seal, high pressure at 

left 


&INPUTS 







TITLE 

= 

' Cylindrical 

step seal. 

high pressure at 

. left' 


I FACE 

= 

0 

ISIUN = 

0 



IGROT 

= 

0 

NOI 

0 

I FLOW 

= 0 

R0 

= 

1 . 0000E+00 

EL 

5.0000E-01 

C 

= 1 . 0000E- 03 

RPM 

= 

5 . 0000E+04 

RPM0 

2 . 5000E+04 

RPMD 

= 0. 0000E+00 

PLEG 

= 

1 . 0000E+03 

PRIG 

0. 0000E+00 

FZD 

= 0. 0000E+00 

vise 

= 

3.0000E-08 

DENS 

1.0000E-04 



EMA 

= 

-2 . 5000E- 01 

ENA 

7 . 9100E- 02 



EMB 

= 

-2 . 5000E-01 

ENB 

7 . 9100E- 02 



HTAP 

= 

0.0000E+00 

HBRL 

0. 0000E + 00 



TOLH 

= 

1.0000E-04 

TOLV 

1.0000E-05 

DUT 

= 1.0000E-06 

I HOME 

= 

0 

NITH 

10 

NITV 

= 30 

NREG 

= 

2 

NRSUB = 

50 50 



ELFR 

= 

5 . 0000E- 01 

5 . 0000E- 01 




ZET 

= 

0.0000E+00 

0 . 0000E+00 




ALP I 

= 

1.0000E+00 

0 . 0000E+00 




BETI 

= 

0 . 0000E+00 

0 . 0000E+00 




DELT 

= 

2 . 0000E-03 

0 . 0000E+00 




NSG 

= 

0 0 





ZETG 

= 

0.0000E+00 

0 . 0000E+00 





/ 


CYLINDRICAL SEAL 

LENGTH, DIAMETER, CLEARANCE = 5.0000E-01, 2.0000E+00, 1.0000E-03 (IN) 

ROTOR, SWIRL AND DIST. SPEEDS = 5.0000E+04, 2.5000E+04, 0.0000E+00 (RPM) 

PRESSURE AT START, END AXIAL BOUNDARIES = 1.0000E+03, 0.0000E+00 (PSI) 

VISCOSITY = 3 . 0000E-08 (PSI-SEC) , DENSITY = 1.0000E-04 (LB-SEC/IN4) 

ERROR CODE = 0, ITERATIONS IN PRIMARY FLOW = 5 

FLOW = 1 . 0041E+01 (IN**3/SEC) 

TORQUE = 5 . 8581E+00 (IN-LB), FILM POWER LOSS = 4.6474E+00 (HP) 

AXIAL REYNOLDS NUMBER = 1.0654E+04 

CIRC. REYNOLDS NUMBERS FOR ROTOR AT SEAL ENDS = 5.2360E+04, 1.7453E+04 


DYNAMIC COEFFICIENTS ( FORCE UNIT / DISP. UNIT ) 


DISP. 

x (IN) 

y (IN) 

phi (RAD) 

psi (RAD) 

FORCE UNIT 

Kx 

1 . 4595E+05 

1. 0148E+05 

- 1 . 4038E+04 

7 . 8450E+04 

LB 

Ky 

- 1 . 0148E+05 

1 . 4595E+05 

- 7 . 8450E+04 

- 1 . 403 8E+04 

LB 

Kphi 

2 . 9812E+03 

-2 . 22 09E+03 

5 . 0418E+03 

9 . 7526E+02 

IN- LB 

Kpsi 

2 . 2209E+03 

2 . 9812E+03 

-9 . 7526E+02 

5 . 0418E+03 

IN-LB 

Bx 

3 . 8776E+01 

9 . 1961E+00 

-1.1539E-01 

5 . 3554E+00 

LB-SEC 

By 

-9 . 1961E+00 

3 . 8776E+01 

-5 . 3554E+00 

-1.1539E-01 

LB-SEC 

Bphi 

1 . 6853E- 01 

-1 . 1378E+00 

3 . 7298E- 01 

5 . 1546E- 02 

IN-LB-SEC 

Bpsi 

1 . 1378E+00 

1 . 6853E- 01 

- 5 . 1546E- 02 

3 . 7298E-01 

IN-LB-SEC 

Ax 

1 . 7549E- 03 

-7 . 6310E- 06 

-3 . 8051E-06 

2 . 2759E-05 

LB-SEC2 

Ay 

7 . 6310E- 06 

1. 7549E-03 

-2 . 2759E- 05 

-3 . 8051E-06 

LB-SEC2 

Aphi 

5 . 1596E- 07 

-3 . 2278E-05 

9 . 7987E- 06 

- 2 . 5418E- 07 

IN-LB-SEC2 

Apsi 

3 . 2278E-05 

5 . 1596E- 07 

2 . 5418E-07 

9.7987E-06 

IN-LB-SEC2 
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( CASE 4 

) 

Cylindrical 

step seal, high pressure at 

right 


&INPUTS 







TITLE 

= 

' Cylindrical 

step seal, 

high pressure at 

: right ' 


I FACE 

= 

0 

ISIUN = 

0 



I GROT 

= 

0 

NOI 

0 

I FLOW = 

0 

R0 

= 

1 . 0000E+00 

EL 

5 . 0000E-01 

C 

1.0000E-03 

RPM 

= 

5 . 0000E+04 

RPM0 

2 . 5000E+04 

RPMD 

0 . 0000E+00 

PLEG 

= 

0.0000E+00 

PRIG 

1. 0000E+03 

FZD 

0 . 0000E+00 

vise 

= 

3.0000E-08 

DENS 

1.0000E-04 



EMA 

= 

-2 . 5000E- 01 

ENA 

7 . 9100E- 02 



EMB 

= 

-2 . 5000E- 01 

ENB 

7 . 9100E- 02 



HTAP 

= 

0 . 0000E+00 

HBRL 

0 . 0000E+00 



TOLH 

= 

1.0000E-04 

TOLV 

1. 0000E-05 

DUT 

1.0000E-06 

I HOME 

= 

0 

NITH 

10 

NITV 

30 

NREG 

= 

2 

NRSUB = 

50 50 



ELFR 

= 

5.0000E-01 

5.0000E-01 




ZET 

= 

0 . 0000E+00 

0 . 0000E+00 




ALPI 

= 

1 . 0000E+00 

0 . 0000E+00 




BETI 

= 

0 . 0000E+00 

0 . 0000E+00 




DELT 

= 

2 . 0000E- 03 

0 . 0000E+00 




NSG 

= 

0 0 





ZETG 

= 

0 . 0000E+00 

0 . 0000E+00 





/ 


CYLINDRICAL SEAL 

LENGTH, DIAMETER, CLEARANCE = 5.0000E-01, 2.0000E+00, 1.0000E-03 (IN) 

ROTOR, SWIRL AND DIST. SPEEDS = 5.0000E+04, 2.5000E+04, 0.0000E+00 (RPM) 

PRESSURE AT START, END AXIAL BOUNDARIES = 0.0000E+00', 1.0000E+03 (PSI) 

VISCOSITY = 3 . 0000E-08 (PSI-SEC) , DENSITY = 1.0000E-04 (LB-SEC/IN4) 

ERROR CODE = 0, ITERATIONS IN PRIMARY FLOW = 3 

FLOW = -1 . 0518E+01 (IN**3/SEC) 

TORQUE = 5 . 8957E+00 (IN-LB), FILM POWER LOSS = 4.6772E+00 (HP) 

AXIAL REYNOLDS NUMBER = 1.1160E+04 

CIRC. REYNOLDS NUMBERS FOR ROTOR AT SEAL ENDS = 5.2360E+04, 1.7453E+04 

DYNAMIC COEFFICIENTS ( FORCE UNIT / DISP. UNIT ) 


DISP. 

x (IN) 

y (IN) 

phi (RAD) 

psi (RAD) 

FORCE UNIT 

Kx 

6 . 6425E+04 

2 . 0961E+04 

7 . 33 02E+03 

-4 . 83 04E+04 

LB 

Ky 

-2 . 0961E+04 

6.6425E+04 

4 . 83 04E+04 

7 . 3302E+03 

LB 

Kphi 

4 . 52 19E+03 

-1 . 8963E+04 

-4.333 5E+03 

3 . 2326E+02 

IN- LB 

Kpsi 

1. 8963E+04 

4 . 5219E+03 

- 3 . 2326E+02 

-4 . 3335E+03 

IN- LB 

Bx 

8 . 5715E+00 

1.387 9E+01 

-5.6371E-01 

-2 . 7546E+00 

LB-SEC 

By 

-1 . 3879E+01 

8 . 5715E+00 

2 . 7 546E+0 0 

-5 . 6371E- 01 

LB -SEC 

Bphi 

4 . 7773E- 01 

- 1 . 7618E+00 

1 . 2625E- 01 

7 . 147 9E- 02 

IN-LB-SEC 

Bpsi 

1 . 7618E+00 

4 . 7773E- 01 

-7 . 1479E-02 

1 . 2625E- 01 

IN-LB-SEC 

Ax 

2 . 5896E-03 

-3 . 1050E- 04 

2 . 4906E- 05 

1 . 0278E-04 

LB-SEC2 

Ay 

3 . 1050E- 04 

2 . 5896E- 03 

-1. 0278E-04 

2 . 4906E- 05 

LB-SEC2 

Aphi 

-1 . 903 0E- 05 

-8 . 7620E- 05 

1 . 336 0E- 05 

-1.5316E-06 

IN-LB-SEC2 

Apsi 

8 . 7 62 0E- 05 

-1. 9030E-05 

1. 5316E-06 

1 . 3360E- 05 

IN-LB-SEC2 
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( CASE 5 

) 

Plain cylindrical seal 

/ 

high pressure 

at left 


&INPUTS 








TITLE 

= 

'Plain cylindrical seal, 

high pressure 

at left' 


I FACE 

= 

0 

ISIUN 

= 

0 



IGROT 

= 

0 

NOI 

= 

0 

I FLOW = 

0 

R0 

= 

1 . 0000E+00 

EL 

= 

5 . 0000E- 01 

C 

1 . 0000E- 03 

RPM 

= 

5 . 0000E+04 

RPM0 

= 

2.5000E+04 

RPMD 

0 . 0000E+00 

PLEG 

= 

1.0000E+03 

PRIG 

= 

0 . 0000E+00 

FZD 

0 . 0000E+00 

vise 

= 

3 . 0000E- 08 

DENS 

= 

1. 0000E-04 



EMA 

= 

-2 . 5000E- 01 

ENA 

= 

7 . 9100E- 02 



EMB 

= 

-2 . 5000E-01 

ENB 

= 

7 . 9100E- 02 



HTAP 

= 

0 . 0000E+00 

HBRL 

= 

0 . 0000E+00 



TOLH 

= 

1.0000E-04 

TOLV 

= 

1. 0000E-05 

DUT 

1.0000E-06 

I HOME 

= 

0 

NITH 

= 

10 

NITV 

30 

NREG 

= 

1 

NRSUB 

= 

100 



ELFR 

= 

1 . 0000E+00 






ZET 

= 

0 . 0000E+00 






ALP I 

= 

0 . 0000E+00 






BETI 

= 

0 . 0000E+00 






DELT 

= 

2 . 0000E- 03 






NSG 

= 

0 






ZETG 

= 

0.0000E+00 







/ 


CYLINDRICAL SEAL 

LENGTH, DIAMETER, CLEARANCE = 5.0000E-01, 2.0000E+00, 1.0000E-03 (IN) 

ROTOR, SWIRL AND DIST. SPEEDS = 5.0000E+04, 2.5000E+04, 0.0000E+00 (RPM) 

PRESSURE AT START, END AXIAL BOUNDARIES = 1.0000E+03, 0.0000E+00 (PSI) 

VISCOSITY = 3.0000E-08 (PSI-SEC) , DENSITY = 1.0000E-04 (LB-SEC/IN4 ) 

ERROR CODE = 0, ITERATIONS IN PRIMARY FLOW = 4 

FLOW = 6 . 2940E+00 (IN**3/SEC) 

TORQUE = 6 . 4988E+00 (IN-LB), FILM POWER LOSS = 5.1557E+00 (HP) 

AXIAL REYNOLDS NUMBER = 6.6781E+03 

CIRC. REYNOLDS NUMBERS FOR ROTOR AT SEAL ENDS = 1.7453E+04, 1.7453E+04 

DYNAMIC COEFFICIENTS ( FORCE UNIT / DISP . UNIT ) 


DISP. 

x (IN) 

y (IN) 

phi (RAD) 

psi (RAD) 

FORCE UNIT 

Kx 

5 . 8071E+04 

2 . 2021E+05 

-1 . 2728E+04 

1 . 63 92E+05 

LB 

Ky 

-2 . 2021E+05 

5 . 8071E+04 

- 1 . 63 92E+05 

-1 . 2728E+04 

LB 

Kphi 

-3 . 7963E+03 

6 . 7330E+03 

- 1 . 8147E+03 

1 . 0783E+03 

IN- LB 

Kpsi 

-6 . 7330E+03 

- 3 . 7963E+03 

- 1 . 0783E+03 

-1 . 8147E+03 

IN- LB 

Bx 

8 . 4163E+01 

1 . 7148E+01 

1 . 5180E- 02 

4 . 86 01E+00 

LB -SEC 

By 

- 1 . 714 8E+01 

8 . 4163E+01 

-4 . 86 01E+00 

1 . 518 0E- 02 

LB-SEC 

Bphi 

-3 . 3856E-02 

1 . 4540E+00 

4 . 1202E- 01 

7 . 2269E- 02 

IN-LB-SEC 

Bpsi 

- 1 . 4540E+00 

-3 . 3856E-02 

-7.2269E-02 

4 . 1202E-01 

IN-LB-SEC 

Ax 

3 . 2700E- 03 

-2 . 7182E- 05 

- 9 . 77 94E- 07 

-2 . 7117E- 06 

LB-SEC2 

Ay 

2 . 7182E- 05 

3 . 2700E- 03 

2 . 7117E- 06 

-9.77 94E- 0 7 

LB-SEC2 

Aphi 

2 . 1509E-06 

6 . 0562E- 06 

1 . 3 7 87E- 05 

-8 . 3503E-08 

IN-LB-SEC2 

Apsi 

-6 . 0562E- 06 

2 . 1509E- 06 

8 . 3503E- 08 

1.378 7E- 05 

IN-LB-SEC2 
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Figure 1 1 Stator with inward pumping grooves for cases 6-8 
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( CASE 6 ) Face seal, grooves pumping inward, high pressure outside 
&INPUTS 


TITLE 

= 

'Face seal, 

grooves pumping inward, 

high’ pressure 

outside ' 

I FACE 

= 

1 

ISIUN = 

0 



I GROT 

= 

0 

NOI 

0 

I FLOW = 

0 

R0 

= 

1 . 0000E+00 

EL 

5 . 0000E- 01 

C 

1.0000E-03 

RPM 

= 

5 . 0000E+04 

RPM0 

2 . 5000E+04 

RPMD 

0 . 0000E+00 

PLEG 

= 

0.0000E+00 

PRIG 

1 . 0000E+03 

FZD 

0 . 0000E+00 

vise 

= 

3 . 0000E-08 

DENS 

1. 0000E-04 



EMA 

= 

-2 . 5000E- 01 

ENA 

7. 9100E-02 



EMB 

= 

- 2 . 5000E- 01 

ENB 

7 . 9100E- 02 



HTAP 

= 

0 . 0000E+00 

HBRL 

0 . 0000E+00 



TOLH 

= 

1.0000E-04 

TOLV 

1. 0000E-05 

DUT 

1. 0000E-06 

I HOME 

= 

0 

NITH 

10 

NITV 

30 

NREG 

= 

2 

NRSUB = 

50 200 



ELFR 

= 

5 . 0000E-01 

5 . 0000E-01 




ZET 

= 

0 . 0000E+00 

0 . 0000E+00 




ALP I 

= 

0 . 0000E+00 

5.0000E-01 




BETI 

= 

0 . 0000E+00 

-2 . 5000E+01 




DELT 

= 

0 . 0000E+00 

2 . 0000E-03 




NSG 

= 

0 28 





ZETG 

= 

0 . 0000E+00 

0 . 0000E+00 





/ 

FACE SEAL 

ID, OD, NOMINAL FILM THICKNESS = 1.0000E+00, 2.0000E+00, 1.0000E-03 (IN) 

ROTOR, SWIRL AND DIST. SPEEDS = 5.0000E+04, 2.5000E+04, 0.0000E+00 (RPM) 

INSIDE, OUTSIDE PRESSURE = 0.0000E+00, 1.0000E+03 (PSI) 

VISCOSITY = 3 . 0000E-08 (PSI-SEC) , DENSITY = 1.0000E-04 (LB-SEC/IN4) 

ERROR CODE = 0, ITERATIONS IN PRIMARY FLOW = 5 

AXIAL LOAD TO BALANCE FACE SEAL = 1.7443E+03 (LB) 


FLOW = - 5 . 5975E+00 (IN**3/SEC) 


TORQUE 

= 2 . 2207E+00 

(IN-LB) , 

FILM POWER : 

LOSS = 1 . 7618E+00 (HP) 

RADIAL 

REYNOLDS NUMBER AT ID, OD 

= 

1.1878E+04, 

5 . 93 91E+03 

CIRC. : 

REYNOLDS NUMBERS FOR ROTOR 

AT ID, OD = 

5 . 8219E+03 , 

3 . 4907E+04 


DYNAMIC COEFFICIENTS ( 

FORCE UNIT / 

DISP. UNIT ) 


DISP. 

z (IN) 

phi (RAD) 

psi (RAD) 

FORCE UNIT 


Kz 

5 . 22 07E+05 

0 . 0000E+00 

0 . 0000E+00 

LB 


Kphi 

0 . 0000E+00 

9 . 0761E+04 

9. 7317E+04 

IN-LB 


Kpsi 

0 . 0000E+00 

- 9 . 7317E+04 

9 . 0761E+04 

IN-LB 


Bz 

1 . 0792E+02 

0 . 0000E+00 

0 . 0000E+00 

LB-SEC • 


Bphi 

0 . 0000E+00 

2 . 73 03E+01 

6 . 5037E+00 

IN-LB-SEC 


Bpsi 

0 . 0000E+00 

-6 . 5037E+00 

2 . 7303E+01 

IN-LB-SEC 


Az 

3 . 5278E-03 

0 . 0000E+00 

0. 0000E+00 

LB-SEC2 


Aphi 

0 . 0000E+00 

9 . 0969E- 04 

- 1 . 3478E- 04 

IN-LB-SEC2 


Apsi 

0 . 0000E+00 

1.3478E-04 

9 . 0969E-04 

IN-LB-SEC2 
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( CASE 7 ) Same as case 6 but with 50% increase in nominal film clearance 


&INPUTS 


TITLE 

= 

' Same as case 6 but with 

50% increase 

in nominal 

film clearance 

I FACE 

= 

1 

ISIUN = 

0 



I GROT 

= 

0 

NOI 

0 

I FLOW = 

0 

R0 

= 

1.0000E+00 

EL 

5 . 0000E-01 

C 

1 . 5000E- 03 

RPM 

= 

5.0000E+04 

RPM0 

2 . 5000E+04 

RPMD 

0 . 0000E+00 

PLEG 

= 

0.0000E+00 

PRIG 

1 . 0000E+03 

FZD 

0 . 0000E+00 

vise 

= 

3 . 0000E-08 

DENS 

1.0000E-04 



EMA 

= 

-2 . 5000E- 01 

ENA 

7 . 9100E- 02 



EMB 

= 

-2 . 5000E-01 

ENB 

7 . 9100E-02 



HTAP 

= 

0 . 0000E+00 

HBRL 

0 . 0000E+00 



TOLH 

= 

1.0000E-04 

TOLV 

1.0000E-05 

DUT 

1.0000E-06 

I HOME 

= 

0 

NITH 

10 

NITV 

30 

NREG 

= 

2 

NRSUB = 

50 200 



ELFR 

= 

5.0000E-01 

5.0000E-01 




ZET 

= 

0 . 0000E+00 

0 . 0000E+00 




ALPI 

= 

0.0000E+00 

5 . 0000E-01 




BETI 

= 

0 . 0000E+00 

-2 . 5000E+01 




DELT 

= 

0 . 0000E+00 

2 . 0000E-03 




NSG 

= 

0 28 





ZETG 

— 

0 . 0000E+00 

0 . 0000E+00 





/ 

FACE SEAL 

ID, OD, NOMINAL FILM THICKNESS = 1.0000E+00, 2.0000E+00, 1.5000E-03 (IN) 

ROTOR, SWIRL AND DIST. SPEEDS = 5.0000E+04, 2.5000E+04, 0 . 0000E+00 (RPM) 

INSIDE, OUTS IDE PRESSURE = 0.0000E+00, 1.0000E+03 (PSI) 

VISCOSITY = 3.0000E-08 (PSI-SEC) , DENSITY = 1.0000E-04 (LB-SEC/IN4) 

ERROR CODE = 0, ITERATIONS IN PRIMARY FLOW = 5 

AXIAL LOAD TO BALANCE FACE SEAL = 1.5888E+03 (LB) 


FLOW = -9 . 1255E+00 (IN**3/SEC) 


TORQUE 

= 1 . 7426E+00 

(IN-LB) , 

FILM POWER : 

LOSS = 1 . 3825E+00 (HP) 

RADIAL 

REYNOLDS NUMBER AT ID , OD 

= 

1 . 9365E+04 , 

9 . 6824E+03 

CIRC. : 

REYNOLDS NUMBERS FOR ROTOR 

AT ID, OD = 

5 . 5391E+03 , 

4 . 3633E+04 


DYNAMIC COEFFICIENTS ( 

FORCE UNIT / 

DISP. UNIT ) 


DISP. 

z (IN) 

phi (RAD) 

psi (RAD) 

FORCE UNIT 


Kz 

1 . 7409E+05 

0 . 0000E+00 

0 . 0000E+00 

LB 


Kphi 

0 . 0000E+00 

1 . 0475E+04 

6 . 1962E+04 

IN-LB 


Kpsi 

0 . 0000E+00 

-6 . 1962E+04 

1. 0475E+04 

IN-LB 


Bz 

6 . 2363E+01 

0 . OOOOE+OO 

0 . 0000E+00 

LB-SEC 


Bphi 

0 . 0000E+00 

1 . 5470E+01 

4 . 4478E+00 

IN-LB-SEC 


Bpsi 

0 . 0000E+00 

-4 . 44 7 8E+00 

1 . 5470E+01 

IN-LB-SEC 


Az 

2 . 3100E- 03 

0. 0000E+00 

0 . 0000E+00 

LB-SEC2 


Aphi 

0 . 0000E+00 

5 . 7624E- 04 

- 7 . 2 022E- 05 

IN-LB-SEC2 


Apsi 

0 . 0000E+00 

7 . 2022E- 05 

5 . 7624E- 04 

IN-LB-SEC2 
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( CASE 8 ) 


Same as case 6 but with IH0ME=1 and 1600 lb. balance load 


& INPUTS 


TITLE 

= 

' Same as case 6 but with 

I HOME =1 and 

1600 lb. 

balance load' 

I FACE 

= 

1 

ISIUN = 

0 




I GROT 

= 

0 

NOI 

0 

I FLOW 

= 

0 

R0 

= 

1 . 0000E+00 

EL 

5 . 0000E-O1 

C 

= 

1.0000E-03 

RPM 

= 

5 . 0000E+04 

RPM0 

2 . 5000E+04 

RPMD 

= 

0 . 0000E+00 

PLEG 

= 

0 . 0000E+00 

PRIG 

1 . 0000E+03 

FZD 

= 

1 . 6000E+03 

vise 

= 

3.0000E-08 

DENS 

1 . 0000E- 04 




EMA 

= 

-2 . 5000E- 01 

ENA 

7 . 9100E- 02 




EMB 

= 

-2 . 5000E- 01 

ENB 

7 . 9100E- 02 




HTAP 

= 

0 . 0000E+00 

HBRL 

0 . 0000E+00 




TOLH 

= 

1 . 0000E- 04 

TOLV 

1.0000E-05 

DUT 

= 

1.0000E-06 

I HOME 

= 

1 

NITH 

10 

NITV 

= 

30 

NREG 

= 

2 

NRSUB = 

50 200 




ELFR 

= 

5.0000E-01 

5.0000E-01 





ZET 

= 

0 . 0000E+00 

0 . 0000E+00 





ALPI 

= 

0 . 0000E+00 

5.0000E-01 





BETI 

= 

0 . 0000E+00 

-2 . 5000E+01 





DELT 

= 

0 . 0000E+00 

2 . 0000E- 03 





NSG 

= 

0 28 






ZETG 

= 

0 . 0000E+00 

0.0000E+00 






/ 

FACE SEAL 

ID, OD, NOMINAL FILM THICKNESS = 1.0000E+00, 2.0000E+00, 1.4389E-03 (IN) 

ROTOR, SWIRL AND DIST. SPEEDS = 5.0000E+04, 2.5000E+04, 0.0000E+00 (RPM) 

INSIDE, OUTSIDE PRESSURE = 0.0000E+00, 1.0000E+03 (PSI) 

VISCOSITY = 3.0000E-08 (PSI-SEC) , DENSITY = 1.0000E-04 (LB-SEC/IN4) 

ERROR CODE = 0, ITERATIONS IN PRIMARY FLOW = 5 

AXIAL LOAD TO BALANCE FACE SEAL = 1.6001E+03 (LB) 

FLOW = -8 . 6909E+00 (IN**3/SEC) 


TORQUE 

= 1 . 7862E+00 

(IN-LB) , 

FILM POWER : 

LOSS = 1 . 4171E+00 (HP) 

RADIAL 

REYNOLDS NUMBER AT ID, OD 

= 

1. 8443E+04, 

9 . 2213E+03 

CIRC. : 

REYNOLDS NUMBERS FOR ROTOR 

AT ID, OD = 

5 . 7012E+03 , 

4 . 2566E+04 


DYNAMIC COEFFICIENTS ( 

FORCE UNIT / 

DISP. UNIT ) 


DISP. 

z (IN) 

phi (RAD) 

psi (RAD) 

FORCE UNIT 


Kz 

1 . 9792E+05 

0. 0000E+00 

0 . 0000E+00 

LB 


Kphi 

0 . 0000E+00 

1 . 5410E+04 

6 . 4918E+04 

IN-LB 


Kpsi 

0 . 0000E+00 

-6 . 4918E+04 

1 . 5410E+04 

IN-LB 


Bz 

6 . 5903E+01 

0 . 0000E+00 

0 . 0000E+00 

LB-SEC 


Bphi 

0 . 0000E+00 

1 . 63 91E+01 

4 . 5984E+00 

IN-LB-SEC 


Bpsi 

0 . 0000E+00 

-4 . 5984E+00 

1 . 6391E+01 

IN-LB-SEC 


Az 

2 . 3 918E- 03 

0. 0000E+00 

0 . 0000E+00 

LB-SEC2 


Aphi 

0 . 0000E+00 

5 . 993 0E- 04 

- 7 . 3444E- 05 

IN-LB-SEC2 


Apsi 

0 . 0000E+00 

7 . 3444E- 05 

5 . 9930E-04 

IN-LB-SEC2 
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( CASE 9 

) 

Childs, Nolan & 

Kilgore Stator 2, 25% 

swirl, 1000 

RPM, NSG=17 

&INPUTS 








TITLE 

= 

' Childs, Nolan 

Sc Kilgore 

Stator 2, 25 

% swirl, 

1000 

RPM, NSG=17 

I FACE 

= 

0 

ISIUN = 

1 




I GROT 

= 

0 

NO I 

0 

I FLOW 

= 

0 

R0 

= 

5 . 0800E- 02 

EL 

5 . 0800E- 02 

C 

= 

3 . 7600E- 04 

RPM 

= 

1. 0000E+03 

RPM0 

2 . 5000E+02 

RPMD 

= 

0 . 0000E+00 

PLEG 

= 

2 . 0000E+06 

PRIG 

0 . 0000E+00 

FZD 

= 

0 . 0000E+00 

vise 

= 

1 . 5400E- 04 

DENS 

1 . 5700E+03 




EMA 

= 

-2 . 5000E- 01 

ENA 

7. 9100E-02 




EMB 

= 

-2 . 5000E- 01 

ENB 

7. 9100E-02 




HTAP 

= 

0 . 0000E+00 

HBRL 

0 . 0000E+00 




TOLH 

= 

1. 0000E-04 

TOLV 

1. 0000E-05 

DUT 

= 

1 . 0000E-06 

I HOME 

= 

0 

NITH 

10 

NITV 

= 

30 

NREG 

= 

1 

NRSUB = 

25 




ELFR 

= 

1.0000E+00 






ZET 

= 

1. 0000E+00 






ALPI 

= 

4 . 9870E-01 






BETI 

= 

1 . 5000E+01 






DELT 

= 

3 . 8000E- 04 






NSG 

= 

17 






ZETG 

= 

0 . 0000E+00 







/ 


CYLINDRICAL SEAL 

LENGTH, DIAMETER, CLEARANCE = 5.0800E-02, 1.0160E-01, 3.7600E-04 (m) 

ROTOR, SWIRL AND DIST. SPEEDS = 1.0000E+03, 2.5000E+02, 0.0000E+00 (RPM) 

PRESSURE AT START, END AXIAL BOUNDARIES = 2.0000E+06, 0 . 0000E+00 (Pa) 

VISCOSITY = 1 . 5400E- 04 (Pa-SEC) , DENSITY = 1.5700E+03 (Kg/m3) 

ERROR CODE = 0, ITERATIONS IN PRIMARY FLOW = 6 

FLOW = 3.3398E-03 (m**3/SEC) 

TORQUE = 2 . 9942E- 02 (N-m) , FILM POWER LOSS = 3.1356E+00 (WATT) 

AXIAL REYNOLDS NUMBER = 2.1335E+05 

CIRC. REYNOLDS NUMBERS FOR ROTOR AT SEAL ENDS = 4.6004E+04, 7.6776E+04 

DYNAMIC COEFFICIENTS ( FORCE UNIT / DISP. UNIT ) 


DISP. 

x (m) 

y (m) 

phi (RAD) 

psi (RAD) 

FORCE UNIT 

Kx 

2 . 6765E+06 

- 3 . 4086E+05 

-1 . 1729E+05 

3 . 9721E+05 

N 

Ky 

3 . 4086E+05 

2 . 6765E+06 

-3 . 9721E+05 

-1.1729E+05 

N 

Kphi 

6 . 1061E+02 

4 . 1684E+04 

-1.2753E+03 

-4 . 4830E+02 

N-m 

Kpsi 

-4 . 1684E+04 

6 . 1061E+02 

4 . 483 0E+02 

-1 . 2753E+03 

N-m 

Bx 

1 . 3 748E+04 

3 . 9289E+03 

5 . 5050E+01 

1 . 6905E+02 

N-SEC 

By 

-3 . 9289E+03 

1 . 3 748E+04 

-1 . 6905E+02 

5 . 5050E+01 

N-SEC 

Bphi 

-1 . 3733E+01 

5 . 9191E+01 

7 . 6324E- 01 

2 . 8576E-01 

N-m-SEC 

Bpsi 

-5 . 9191E+01 

- 1 . 3 73 3E+01 

-2 . 8576E- 01 

7 . 6324E- 01 

N-m-SEC 

Ax 

5 . 5942E+00 

- 1 . 3283E+00 

-2 . 8954E-02 

-3 . 1742E- 02 

N-SEC2 

Ay 

1 . 3283E+00 

5 . 5942E+00 

3 . 1742E- 02 

-2 . 8954E- 02 

N-SEC2 

Aphi 

4 . 7119E- 03 

7 . 43 99E- 03 

3 . 3120E-04 

-6 . 973 0E- 05 

N-m-SEC2 

Apsi 

-7 . 4399E-03 

4 . 7119E- 03 

6 . 973 0E- 05 

3 . 3120E-04 

N-m-SEC2 
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( CASE 10 ) 


Same as Case 9 except NSG=0 


& INPUTS 


TITLE 

= 

' Same as Case 

9 except NSG=0' 



I FACE 

= 

0 

ISIUN = 

1 



I GROT 

= 

0 

NOI 

0 

I FLOW 

= 0 

R0 

= 

5 . 0800E- 02 

EL 

5 . 0800E- 02 

C 

= 3 . 7600E- 04 

RPM 

= 

1.0000E+03 

RPM0 

2 . 5000E+02 

RPMD 

= 0.0000E+00 

PLEG 

= 

2 . 0000E+06 

PRIG 

0.0000E+00 

FZD 

= 0.0000E+00 

vise 

= 

1 . 5400E- 04 

DENS 

1 . 5700E+03 



EMA 

= 

-2 . 5000E- 01 

ENA 

7 . 9100E- 02 



EMB 

= 

- 2 . 5000E- 01 

ENB 

7 . 9100E- 02 



HTAP 

= 

0.0000E+00 

HBRL 

0.0000E+00 



TOLH 

= 

1.0000E-04 

TOLV 

1. 0000E-05 

DUT 

= 1.0000E-06 

I HOME 

= 

0 

NITH 

10 

NITV 

= 30 

NREG 

= 

1 

NRSUB = 

25 



ELFR 

= 

1 . 0000E+00 





ZET 

= 

1. 0000E+00 





ALP I 

= 

4 . 9870E- 01 





BETI 

= 

1 . 5000E+01 





DELT 

= 

3 . 8000E- 04 





NSG 

= 

0 





ZETG 

— 

0.0000E+00 






/ 

CYLINDRICAL SEAL 

LENGTH, DIAMETER, CLEARANCE = 5.0800E-02, 1.0160E-01, 3.7600E-04 (m) 

ROTOR, SWIRL AND DIST. SPEEDS = 1.0000E+03, 2.5000E+02, O.OOOOE+OO (RPM) 

PRESSURE AT START, END AXIAL BOUNDARIES = 2.0000E+06, 0.0000E+00 (Pa) 

VISCOSITY = 1 . 5400E-04 (Pa-SEC) , DENSITY = 1.5700E+03 (Kg/m3) 

ERROR CODE = 0, ITERATIONS IN PRIMARY FLOW = 5 

FLOW = 5.2060E-03 (m**3/SEC) 

TORQUE = 4 . 3118E-01 (N-m) , FILM POWER LOSS = 4 . 5153E+01 (WATT) 

AXIAL REYNOLDS NUMBER = 3.3256E+05 

CIRC. REYNOLDS NUMBERS FOR ROTOR AT SEAL ENDS = 4.6004E+04, 1.9099E+04 

DYNAMIC COEFFICIENTS ( FORCE UNIT / DISP . UNIT ) 


DISP. 

x (m) 

y (m) 

phi (RAD) 

psi (RAD) 

FORCE UNIT 

Kx 

5 . 7006E+06 

-4 . 4186E+05 

-3 . 8223E+04 

4 . 7259E+05 

N 

Ky 

4 . 4186E+05 

5 . 7006E+06 

-4 . 7259E+05 

-3 . 8223E+04 

N 

Kphi 

1 . 5121E+03 

5 . 4046E+04 

-2 . 9642E+03 

-1. 8965E+02 

N-m 

Kpsi 

- 5 . 4046E+04 

1 . 512 1E+03 

1 . 8965E+02 

-2 . 9642E+03 

N-m 

Bx 

1 . 1165E+04 

1.228 9E+03 

1 . 6233E+01 

2 . 4667E+02 

N-SEC 

By 

- 1 . 2289E+03 

1. 1165E+04 

-2 . 4667E+02 

1 . 6233E+01 

N-SEC 

Bphi 

-6 . 2 97 9E+00 

8 . 1492E+01 

1 . 9357E- 01 

1 . 5111E- 01 

N-m-SEC 

Bpsi 

-8 . 14 92E+01 

-6 . 2979E+00 

- 1 . 5111E- 01 

1. 9357E-01 

N-m-SEC 

Ax 

6 . 3347E+00 

- 3 . 42 09E- 01 

-7 . 1525E- 03 

-4 . 4958E- 02 

N-SEC2 

Ay 

3 . 4209E- 01 

6 . 3347E+00 

4 . 4 958E- 02 

-7 . 1525E-03 

N-SEC2 

Aphi 

2 . 6813E- 03 

1.7398E-02 

6 . 0569E- 04 

-5.3545E-05 

N-m-SEC2 

Apsi 

-1. 7398E-02 

2.6813E-03 

5 . 3545E-05 

6 . 0569E-04 

N-m- SEC2 
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5. Verification 


An available MTI computer code has been used to obtain a partial check on the results of SPIRALI. This 
code implements a two dimensional finite difference solution to the spiral groove equations that is quite 
different from the perturbation methods used in SPIRALI. It does not treat inertia effects and uses the Ng 
and Pan (ref. 12) linearized turbulence theory which limits its validity to cases where turbulence is 
dominated by circumferential component of the flow. The assumptions used in the MTI code will be 
identical to those used in SPIRALI when the flow is laminar and inertia effects are neglected. The bulk 
flow model used in SPIRALI should also produce similar (although not identical) results to those obtained 
with the MTI code under turbulent conditions when inertia effects are neglected, there is no significant 
seal pressure difference, the circumferential Reynolds numbers are between 10* and 1GP and the 
transverse Reynolds numbers are small compared with the circumferential ones. 

Four test cases obtained with SPIRALI are presented for comparison with the MTI code on the following 
pages. Cases 1 and 3 correspond to laminar flow in a spiral grooved cylindrical and face seal 
respectively. It can be seen that the results of both cases agree extremely well with the results of the 
MTI code. The small discrepancies can easily be explained by truncation error associated with the finite 
difference solution used in the MTI code. Cases 2 and 4 provide comparisons similar to cases 1 and 3 
under conditions of turbulent flow. It can be seen that the Reynolds numbers are in the general range of 
validity of the MTI code and all results agree to within 10% for both cases which is acceptable 
considering the differences in the turbulent flow models. 

While the above cases provide a partial check on turbulent, spiral groove seals, they do not provide any 
validation on the treatment of inertia effects or seal behavior under high imposed pressure gradients. In 
order to obtain verification of the treatment of these effects, comparisons are presented between the 
results published by Childs in reference 4 and those obtained with SPIRALI for a plain ungrooved 
cylindrical seal. The output from SPIRALI is given as cases 5 - 8 on the following pages and 
comparisons of flow rates and rotordynamic coefficients with those published by Childs are given in 
Table 1. Although agreement is fairly good, discrepancies in results of up to 10%, such as that which 
occurs in the direct stiffness for case 6 where 4 n =r 0 u/2 and L/D=1, are somewhat surprising since the 
assumptions for that case are identical with those used by Childs. In order to obtain a further 
independent check on the results, a reduced form of the differential equations published by Childs for 
predicting the direct and cross stiffness coefficients when 4 n =r 0 u)/ 2 , were programmed and solved by the 
Runge-Kutta method. The program produced values of K xx , and Q that agree to four places with the 
case 6 results produced by SPIRALI. It is thus likely that the small discrepancies between results 
obtained from SPIRALI and those published in reference 4 are due to truncation errors in the latter 
solution. Additional verification for helical groove seals is included in Appendix A. 
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( CASE 1 

) 

Cylindrical 

seal with grooves, laminar, no press. 

grad. 

&INPUTS 







TITLE 

= 

' Cylindrical 

seal with grooves, laminar, no press, 

. grad. ' 

I FACE 

= 

0 

ISIUN = 

0 



I GROT 

= 

0 

NOI 

2 I FLOW 

= 

1 

R0 

= 

1 . 0000E+00 

EL 

5 . 0000E- 01 C 

= 

1.0000E-03 

RPM 

= 

5 . 0000E+04 

RPM0 

2 . 5000E+04 RPMD 

= 

0.0000E+00 

PLEG 

= 

0 . 0000E+00 

PRIG 

0.0000E+00 FZD 

= 

0 . 0000E+00 

vise 

= 

3.0000E-08 

DENS 

0.0000E+00 



EMA 

= 

-2 . 5000E-01 

ENA 

7 . 9100E- 02 



EMB 

= 

- 2 . 5000E- 01 

ENB 

7 . 9100E- 02 



HTAP 

= 

0.0000E+00 

HBRL 

0 . 0000E+00 



TOLH 

= 

1.0000E-04 

TOLV 

1.0000E-05 DUT 

= 

1.0000E-06 

I HOME 

= 

0 

NITH 

10 NITV 

= 

30 

NREG 

= 

2 

NRSUB = 

50 50 



ELFR 

= 

5.0000E-01 

5.0000E-01 




ZET 

= 

0.0000E+00 

0 . 0000E+00 




ALP I 

= 

5.0000E-01 

0.0000E+00 




BETI 

= 

2 . 5000E+01 

0.0000E+00 




DELT 

= 

2.0000E-03 

0.0000E+00 




NSG 

= 

0 0 





ZETG 

= 

0.0000E+00 

0.0000E+00 





/ 


CYLINDRICAL SEAL, INERTIA NEGLECTED 

LENGTH, DIAMETER, CLEARANCE = 5.0000E-01, 2.0000E+00, 1.0000E-03 (IN) 

ROTOR, SWIRL AND DIST. SPEEDS = 5.0000E+04, 2.5000E+04, O.OOOOE+OO (RPM) 

PRESSURE AT START, END AXIAL BOUNDARIES = 0.0000E+00, 0.0000E+00 (PSI) 

VISCOSITY = 3 . 0000E-08 (PSI-SEC) , DENSITY = 0.0000E+00 (LB-SEC/IN4) 

ERROR CODE = 0, ITERATIONS IN PRIMARY FLOW = 2 

FLOW = 1 . 1506E+00 (IN**3/SEC) 

TORQUE = 4 . 3909E-01 (IN-LB) , FILM POWER LOSS = 3.4835E-01 (HP) 

AXIAL REYNOLDS NUMBER = 0.0000E+00 

CIRC. REYNOLDS NUMBERS FOR ROTOR AT SEAL ENDS = 0.0000E+00, 0.0000E+00 

DYNAMIC COEFFICIENTS ( FORCE UNIT / DISP . UNIT ) 


DISP. 

x (IN) 

y (IN) 

phi (RAD) 

psi (RAD) 

FORCE UNIT 

Kx 

2 . 7021E+04 

1 . 3511E+04 

-6 . 3081E+02 

-1 . 2959E+03 

LB 

Ky 

-1 . 3511E+04 

2 . 7021E+04 

1 . 2 959E+03 

-6 . 3081E+02 

LB 

Kphi 

2 . 3387E+02 

-6 . 5242E+01 

1 . 66 99E+02 

8 . 6419E+01 

IN- LB 

Kpsi 

6 . 5242E+01 

2 . 3387E+02 

- 8 . 6419E+01 

1 . 6699E+02 

IN- LB 

Bx 

5 . 1297E+00 

7 . 8477E- 11 

-2 . 0242E-02 

1 . 3 973E- 01 

LB -SEC 

By 

-7 . 8477E- 11 

5 . 1297E+00 

- 1 . 3 973E- 01 

-2 . 0242E- 02 

LB -SEC 

Bphi 

-2 . 0242E- 02 

- 1 . 3 973E- 01 

3 . 0437E-02 

3 . 8319E- 13 

IN-LB-SEC 

Bpsi 

1 . 3973E-01 

-2 . 0242E- 02 

-3 . 8319E- 13 

3 . 0437E-02 

IN-LB-SEC 

Ax 

4 . 7708E- 19 

4 . 0552E- 18 

7 . 4544E-20 

-1.7891E-19 

LB-SEC2 

Ay 

-4 . 0552E- 18 

4 . 7708E- 19 

1 . 7891E-19 

7 . 4544E-20 

LB-SEC2 

Aphi 

1 . 1927E- 19 

5 . 4045E- 2 0 

- 1 . 8636E-20 

3 . 3545E- 20 

IN-LB-SEC2 

Apsi 

- 5 . 4045E-20 

1 . 1927E- 19 

-3 . 3545E-20 

- 1 . 8636E- 2 0 

IN-LB-SEC2 
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RESULTS OF MTI SPIRAL GROOVE CODE FOR COMPARISON WITH CASE 1 


FLOW = 1 . 1506E+00 (IN**3/SEC) 

TORQUE = 4.3909E-01 (IN-LB), FILM POWER LOSS = 3.4835E-01 (HP) 

DYNAMIC COEFFICIENTS ( FORCE UNIT / DISP. UNIT ) 


DISP. 

x (IN) 

y (IN) 

phi (RAD) 

psi (RAD) 

FORCE UNIT 

Kx 

2 .6996E+04 

1 . 3497E+04 

-6.2966E+02 

-1 . 2947E+03 

LB 

Ky 

- 1 . 3497E+04 

2 . 6996E+04 

1 . 2947E+03 

-6 . 2966E+02 

LB 

Kphi 

2 . 3429E+02 

-6 . 4947E+01 

1 . 6666E+02 

8 . 6250E+01 

IN-LB 

Kpsi 

6 . 4947E+01 

2 . 342 9E+02 

- 8 . 6250E+01 

1 . 6666E+02 

IN- LB 

Bx 

5 . 1280E+00 

2 . 4507E- 07 

-2 . 0178E- 02 

1 . 3967E- 01 

LB-SEC 

By 

9 . 3415E- 09 

5.1280E+00 

-1.3967E-01 

-2 . 0178E- 02 

LB-SEC 

Bphi 

-2 . 0178E- 02 

-1 . 3967E-01 

3 . 0411E- 02 

3 . 0462E- 09 

IN-LB-SEC 

Bpsi 

1.3967E-01 

-2 . 0178E- 02 

4 . 9311E- 09 

3 . 0411E- 02 

IN-LB-SEC 
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( CASE 2 

) 

Cylindrical 

seal with grooves, turbulent 

. , no press, grad. 

& INPUTS 







TITLE 

= 

* Cylindrical 

seal with grooves, turbulent, no 

press, grad.' 

I FACE 

= 

0 

ISIUN = 

0 



I GROT 

= 

0 

NO I 

2 

I FLOW 

= 1 

R0 

= 

1 . 0000E+00 

EL 

5.0000E-01 

C 

= 1.0000E-03 

RPM 

= 

5.0000E+04 

RPMO 

2 . 5000E+04 

RPMD 

= 0.0000E+00 

PLEG 

= 

0.0000E+00 

PRIG 

0.0000E+00 

FZD 

= 0.0000E+00 

vise 

= 

3 . 0000E-08 

DENS 

1.0000E-04 



EMA 

= 

-2 . 5000E- 01 

ENA 

7 . 9100E- 02 



EMB 

= 

-2 . 5000E- 01 

ENB 

7 . 9100E- 02 



HTAP 

= 

0.0000E+00 

HBRL 

0.0000E+00 



TOLH 

= 

1.0000E-04 

TOLV 

1.0000E-05 

DUT 

= 1.0000E-06 

I HOME 

= 

0 

NITH 

10 

NITV 

= 30 

NREG 

= 

2 

NRSUB = 

50 50 



ELFR 

= 

5.0000E-01 

5.0000E-01 




ZET 

= 

0.0000E+00 

0.0000E+00 




ALPI 

= 

5.0000E-01 

0 . 0000E+00 




BETI 

= 

2 . 5000E+01 

0.0000E+00 




DELT 

= 

2 . 0000E-03 

0.0000E+00 




NSG 

= 

0 0 





ZETG 

= 

0.0000E+00 

0.0000E+00 





/ 


CYLINDRICAL SEAL, INERTIA NEGLECTED 

LENGTH, DIAMETER, CLEARANCE = 5.0000E-01, 2.0000E+00, 1.0000E-03 (IN) 

ROTOR, SWIRL AND DIST. SPEEDS = 5.0000E+04, 2.5000E+04, 0.0000E+00 (RPM) 

PRESSURE AT START, END AXIAL BOUNDARIES = 0.0000E+00, O.OOOOE+OO (PSI) 

VISCOSITY = 3 . 0000E-08 (PSI-SEC) , DENSITY = 1.0000E-04 (LB-SEC/IN4) 

ERROR CODE = 0, ITERATIONS IN PRIMARY FLOW = 3 

FLOW = 1 . 7348E+00 (IN**3/SEC) 

TORQUE = 6 . 1300E+00 (IN-LB), FILM POWER LOSS = 4.8631E+00 (HP) 

AXIAL REYNOLDS NUMBER = 1.8407E+03 

CIRC. REYNOLDS NUMBERS FOR ROTOR AT SEAL ENDS = 3.8561E+04, 1.7453E+04 

DYNAMIC COEFFICIENTS ( FORCE UNIT / DISP. UNIT ) 


DISP. 

x (IN) 

y (IN) 

phi (RAD) 

psi (RAD) 

FORCE UNIT 

Kx 

2 . 0689E+05 

1. 0073E+05 

-1. 9681E+03 

- 9 . 2 019E+03 

LB 

Ky 

- 1 . 0073E+05 

2 . 0689E+05 

9 . 2 019E+03 

- 1 . 9681E+03 

LB 

Kphi 

8 . 8445E+02 

-2 . 4291E+02 

1 . 1123E+03 

4 . 9214E+02 

IN- LB 

Kpsi 

2 . 42 91E+02 

8 . 8445E+02 

-4 . 9214E+02 

1 . 1123E+03 

IN- LB 

Bx 

3 . 7161E+01 

3 . 3 996E- 07 

-6 . 9734E- 02 

5 . 2 914E- 01 

LB-SEC 

By 

-3 . 3996E-07 

3 . 7161E+01 

- 5 . 2 914E- 01 

-6 . 9734E- 02 

LB-SEC 

Bphi 

-6 . 9734E- 02 

- 5 . 2 914E- 01 

1. 7825E-01 

2 . 1459E- 12 

IN-LB-SEC 

Bpsi 

5 . 2914E-01 

-6 . 9734E- 02 

-2 . 1459E- 12 

1. 7825E-01 

IN-LB-SEC 

Ax 

3 . 8167E-18 

7 . 6333E- 18 

2 . 3854E-18 

-1 . 6698E- 18 

LB-SEC2 

Ay 

-7 . 6333E-18 

3 . 8167E-18 

1 . 6698E-18 

2 . 3854E-18 

LB-SEC2 

Aphi 

1 . 2 076E- 18 

-3 . 5781E- 19 

1 . 4909E-19 

-4 . 4727E-20 

IN-LB-SEC2 

Apsi 

3 . 5781E-19 

1 . 2076E-18 

4 . 4727E-20 

1 . 4909E- 19 

IN-LB-SEC2 
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RESULTS OF MTI SPIRAL GROOVE CODE FOR COMPARISON WITH CASE 2 


FLOW = 

1 . 7076E+00 

(IN**3/ SEC) 




TORQUE 

= 5 . 8 953E+00 

(IN-LB) , 

FILM POWER 

LOSS = 4 . 6770E+00 (HP) 


DYNAMIC COEFFICIENTS ( 

FORCE UNIT / DISP. UNIT ) 


DISP. 

x (IN) 

y (in) 

phi (RAD) 

psi (RAD) 

FORCE UNIT 

Kx 

2 . 1565E+05 

1 . 0203E+05 

-2 . 0751E+03 

- 9 . 863 7E+03 

LB 

Ky 

-1 . 0203E+05 

2 . 156 5E+05 

9 . 8637E + 03 

-2 . 0 7 51E + 03 

LB 

Kphi 

9 . 7231E+02 

-2 . 4006E+02 

1 . 1627E+03 

5 . 0513E+02 

IN -LB 

Kpsi 

2 . 4006E+02 

9 . 7231E+02 

- 5 . 0513E+02 

1 . 1627E+03 

IN-LB 

Bx 

3 . 7602E+01 

2 . 8881E- 06 

-6.8 925E- 02 

5.6836E-01 

LB -SEC 

By 

3 . 6843E- 07 

3 . 7602E+01 

- 5 . 6836E- 01 

-6 . 8930E-02 

LB -SEC 

Bphi 

-6 . 8929E- 02 

- 5 . 6836E- 01 

1.8265E-01 

6 . 4583E- 08 

IN-LB-SEC 

Bpsi 

5 . 6836E- 01 

-6 . 8929E- 02 

-3 . 9709E- 09 

1. 8265E-01 

IN-LB-SEC 
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( CASE 3 ) 


Face seal with grooves, laminar, no press, grad. 


& INPUTS 


TITLE 

= 

' Face seal 

with grooves, 

laminar, no 

press, grad. 

i 

I FACE 

= 

1 

ISIUN = 

0 




I GROT 

= 

0 

NOI 

2 

I FLOW 

= 

1 

R0 

= 

1 . 0000E+00 

EL 

5.0000E-01 

C 

= 

1.0000E-03 

RPM 

= 

5 . 0000E+04 

RPM0 

2 . 5000E+04 

RPMD 

= 

0 . 0000E+00 

PLEG 

= 

0.0000E+00 

PRIG 

0 . 0000E+00 

FZD 

= 

0 . 0000E+00 

vise 

= 

3 . 0000E-08 

DENS 

0 . 0000E+00 




EMA 

= 

-2 . 5000E-01 

ENA 

7 . 9100E- 02 




EMB 

= 

-2 . 5000E- 01 

ENB 

7 . 9100E- 02 




HTAP 

= 

0 . 0000E+00 

HBRL 

0 . 0000E+00 




TOLH 

= 

1. 0000E-04 

TOLV 

1. 0000E-05 

DUT 

= 

1.0000E-06 

I HOME 

= 

0 

NITH 

10 

NITV 

= 

30 

NREG 

= 

2 

NRSUB = 

50 50 




ELFR 

= 

5 . 0000E- 01 

5.0000E-01 





ZET 

= 

0 . 0000E+00 

0 . 0000E+00 





ALP I 

= 

5 . 0000E-01 

0 . 0000E+00 





BETI 

= 

2 . 5000E+01 

0 . 0000E+00 





DELT 

= 

2 . 0000E-03 

0 . 0000E+00 





NSG 

= 

0 0 






ZETG 

= 

0 . 0000E+00 

0 . 0000E+00 






/ 

FACE SEAL, INERTIA NEGLECTED 

ID, OD, NOMINAL FILM THICKNESS = 1.0000E+00, 2.0000E+00, 1.0000E-03 (IN) 

ROTOR, SWIRL AND DIST. SPEEDS = 5.0000E+04, 2.5000E+04, 0.0000E+00 (RPM) 

INSIDE, OUTSIDE PRESSURE = 0.0000E+00, 0 . 0000E+00 (PSI) 

VISCOSITY = 3 . 0000E-08 (PSI-SEC) , DENSITY = 0.0000E+00 (LB-SEC/IN4) 

ERROR CODE = 0, ITERATIONS IN PRIMARY FLOW = 2 

AXIAL LOAD TO BALANCE FACE SEAL = 1.0391E+01 (LB) 

FLOW = 5 . 7837E-01 (IN**3/SEC) 


TORQUE 

= 2 . 1724E- 01 

(IN-LB) , 

FILM POWER : 

LOSS = 1 . 7234E- 01 (HP) 

RADIAL 

REYNOLDS NUMBER AT ID, OD 

— 

0 . 0000E+00 , 0 . 0000E+00 

CIRC. : 

REYNOLDS NUMBERS FOR ROTOR 
DYNAMIC COEFFICIENTS ( 

AT ID, OD = 
FORCE UNIT / 

0 . 0000E+00 , 0 . 0000E+00 

DISP. UNIT ) 

DISP. 

z (IN) 

phi (RAD) 

psi (RAD) 

FORCE UNIT 

Kz 

2 . 3640E+04 

0 . 0000E+00 

0 . 0000E+00 

LB 

Kphi 

0 . 0000E+00 

6 . 2646E+03 

7 . 0922E+03 

IN- LB 

Kpsi 

0 . 0000E+00 

-7 . 0922E+03 

6.2646E+03 

IN-LB 

Bz 

8 . 9011E+00 

0 . 0000E+00 

0 . 0000E+00 

LB-SEC 

Bphi 

0 . 0000E+00 

2 . 5921E+00 

2 . 6323E- 05 

IN-LB-SEC 

Bpsi 

0 . 0000E+00 

-2 . 6323E- 05 

2 . 5921E+00 

IN-LB-SEC 

Az 

0 . 0000E+00 

0 . 0000E+00 

0 . 0000E+00 

LB-SEC2 

Aphi 

0 . 0000E+00 

-1.1927E-19 

1 . 1927E- 19 

IN-LB-SEC2 

Apsi 

0 . 0000E+00 

-1.1927E-19 

- 1 . 1927E- 19 

IN-LB-SEC2 
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RESULTS OF MTI SPIRAL GROOVE CODE FOR COMPARISON WITH CASE 3 


AXIAL LOAD TO BALANCE FACE SEAL = 1.0391E+01 (LB) 

FLOW = 5 . 7837E-01 (IN**3/SEC) 

TORQUE = 2 . 1724E-01 (IN-LB), FILM POWER LOSS = 1.7234E-01 


DYNAMIC COEFFICIENTS ( FORCE UNIT / DISP. UNIT ) 


DISP. 

z (IN) 

phi (RAD) 

psi (RAD) 

FORCE UNIT 

Kz 

2 . 3639E+04 

7 . 1718E- 03 

7 . 1755E- 03 

LB 

Kphi 

5.8063E-06 

6 . 2581E+03 

7 . 0834E+03 

IN-LB 

Kpsi 

1 . 1824E- 05 

-7 . 0834E+03 

6 . 2581E+03 

IN- LB 

Bz 

8 . 9022E+00 

-2 . 7999E- 07 

-3 . 9768E- 07 

LB-SEC 

Bphi 

-6 . 6333E- 08 

2 . 5909E+00 

2 . 9199E- 06 

IN-LB-SEC 

Bpsi 

- 1 . 9181E- 08 

-2 . 9131E- 06 

2 . 5909E+00 

IN-LB-SEC 


(HP) 
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( CASE 4 ) 


Face seal with grooves, turbulent, no press, grad. 


& INPUTS 


TITLE 

= 

1 Face seal 

with grooves, 

turbulent, no 

press . 

grad. * 

I FACE 

= 

1 

ISIUN = 

0 




I GROT 

= 

0 

NOI 

2 

I FLOW 

= 

1 

R0 

= 

1. 0000E+00 

EL 

5.0000E-01 

C 

= 

1.0000E-03 

RPM 

= 

5.0000E+04 

RPM0 

2 . 5000E+04 

RPMD 

= 

0 . 0000E+00 

PLEG 

= 

0.0000E+00 

PRIG 

O.OOOOE+OO 

FZD 

= 

0.0000E+00 

vise 

= 

3 . 0000E- 08 

DENS 

1.0000E-04 




EMA 

= 

-2 . 5000E- 01 

ENA 

7 . 9100E- 02 




EMB 

= 

- 2 . 5000E- 01 

ENB 

7 . 9100E- 02 




HTAP 

= 

0.0000E+00 

HBRL 

O.OOOOE+OO 




TOLH 

= 

1.0000E-04 

TOLV 

1.0000E-05 

DUT 

= 

1.0000E-06 

I HOME 

= 

0 

NITH 

10 

NITV 

= 

30 

NREG 

= 

2 

NRSUB = 

50 50 




ELFR 

= 

5.0000E-01 

5 . 0000E-01 





ZET 

= 

0.0000E+00 

0.0000E+00 





ALP I 

= 

5.0000E-01 

0.0000E+00 





BETI 

= 

2 . 5000E+01 

0.0000E+00 





DELT 

= 

2 . 0000E-03 

0.0000E+00 





NSG 

= 

0 0 






ZETG 

= 

0.0000E+00 

0.0000E+00 






FACE SEAL, INERTIA NEGLECTED 

ID, OD, NOMINAL FILM THICKNESS = 1.0000E+00, 2.0000E+00, 1.0000E-03 (IN) 

ROTOR, SWIRL AND DIST. SPEEDS = 5.0000E+04, 2.5000E+04, O.OOOOE+OO (RPM) 

INSIDE, OUTSIDE PRESSURE = 0.0000E+00, 0.0000E+00 (PSI) 

VISCOSITY = 3 . 0000E-08 (PSI-SEC) , DENSITY = 1.0000E-04 (LB-SEC/IN4) 

ERROR CODE = 0, ITERATIONS IN PRIMARY FLOW = 3 

AXIAL LOAD TO BALANCE FACE SEAL = 5 . 8223E+01 (LB) 


FLOW = 

7 . 2459E- 01 

(IN**3/ SEC) 




TORQUE 

= 2 . 4948E+00 

(IN-LB) , 

FILM POWER : 

LOSS = 1 . 9792E+00 (HP) 

RADIAL 

REYNOLDS NUMBER AT ID, OD 

= 

1 . 5376E+03 , 

7 . 6881E+02 

CIRC. REYNOLDS NUMBERS FOR ROTOR 

AT ID, OD = 

1 . 8 985E+04 , 

1 . 7453E+04 


DYNAMIC COEFFICIENTS ( 

FORCE UNIT / 

DISP. UNIT ) 


DISP. 

z (IN) 

phi (RAD) 

psi (RAD) 

FORCE UNIT 


Kz 

1 . 2 74 0E+05 

0.0000E+00 

0 . 0000E+00 

LB 


Kphi 

0 . 0000E+00 

3 . 4989E+04 

3 . 86 92E+04 

IN-LB 


Kpsi 

0 . 0000E+00 

-3 . 8692E+04 

3 . 4 98 9E + 04 

IN-LB 


Bz 

4 . 8105E+01 

0 . 0000E+00 

0 . 0000E+00 

LB-SEC 


Bphi 

0 . 0000E+00 

1 . 4225E+01 

6 . 4594E- 05 

IN-LB-SEC 


Bpsi 

0 . 0000E+00 

-6 .4594E-05 

1 . 4225E+01 

IN-LB-SEC 


Az 

0 . 0000E+00 

0.0000E+00 

0 . 0000E+00 

LB-SEC2 


Aphi 

0 . 0000E+00 

1. 9083E-18 

- 9 . 5417E- 19 

IN-LB-SEC2 


Apsi 

0 . 0000E+00 

9 . 5417E- 19 

1. 9083E-18 

IN-LB-SEC2 
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RESULTS OF MTI SPIRAL GROOVE CODE FOR COMPARISON WITH CASE 4 
AXIAL LOAD TO BALANCE FACE SEAL = 5.8223E+01 (LB) 

FLOW = 7.2032E-01 (IN**3/SEC) 

TORQUE = 2 . 4013E+00 (IN-LB), FILM POWER LOSS = 1.9050E+00 

DYNAMIC COEFFICIENTS ( FORCE UNIT / DISP. UNIT ) 

DISP. Z (IN) phi (RAD) psi (RAD) FORCE UNIT 

Kz 1 . 3685E+05 5.1323E-02 5.1174E-02 LB 

Kphi 8 . 1427E-06 3.7351E+04 4.0251E+04 IN-LB 

Kpsi 3 . 8131E-05 -4 . 0251E+04 3.7351E+04 IN-LB 

Bz 4 . 9999E+01 -3.3078E-06 -3.4927E-06 LB-SEC 

Bphi 9 . 6839E-07 1.4768E+01 7.7653E-06 IN-LB-SEC 

Bpsi 2 . 8875E-07 -6.9217E-06 1.4768E+01 IN-LB-SEC 


(HP) 
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( CASE 5 

) 

Childs finite 

length solution, RPM0= 

RPM/ 2 , L/D= . 2 


&INPUTS 








TITLE 

= 

'Childs finite 

length 

solution, RPM0 

=RPM/2, L/D=. 

2 ' 

I FACE 

= 

0 

ISIUN 

= 

1 



I GROT 

= 

0 

NO I 

= 

0 

I FLOW = 

0 

R0 

= 

7.6200E-02 

EL 

= 

3 . 048 0E- 02 

C 

1 . 9050E- 04 

RPM 

= 

3 . 6000E+03 

RPM0 

= 

1 . 8000E+03 

RPMD 

0 . 0000E+00 

PLEG 

= 

3 . 4400E+06 

PRIG 

= 

0 . 0000E+00 

FZD 

0 . 0000E+00 

vise 

= 

1 . 2950E-03 

DENS 

= 

1. 0000E+03 



EMA 

= 

-2 . 5000E-01 

ENA 

= 

7 . 9100E- 02 



EMB 

= 

-2 . 5000E- 01 

ENB 

= 

7 . 9100E- 02 



HTAP 

= 

0 . 0000E+00 

HBRL 

= 

0 . 0000E+00 



TOLH 

= 

1. 0000E-04 

TOLV 

= 

1. 0000E-05 

DUT 

1 . 0000E-06 

I HOME 

= 

0 

NITH 

= 

10 

NITV 

30 

NREG 

= 

1 

NRSUB 

= 

200 



ELFR 

= 

1. 0000E+00 






ZET 

= 

1 . 0000E- 01 






ALP I 

= 

0.0000E+00 






BETI 

= 

0 . 0000E+00 






DELT 

= 

0 . 0000E+00 






NSG 

= 

0 






ZETG 

= 

0.0000E+00 







/ 


CYLINDRICAL SEAL 

LENGTH, DIAMETER, CLEARANCE = 3.0480E-02, 1.5240E-01, 1.9050E-04 (m) 

ROTOR, SWIRL AND DIST. SPEEDS = 3.6000E+03, 1.8000E+03, 0.0000E+00 (RPM) 

PRESSURE AT START, END AXIAL BOUNDARIES = 3.4400E+06, 0.0000E+00 (Pa) 

VISCOSITY = 1 . 2950E-03 (Pa-SEC) , DENSITY = 1.0000E+03 (Kg/m3) 

ERROR CODE = 0, ITERATIONS IN PRIMARY FLOW = 4 

FLOW = 4.0061E-03 (m**3/SEC) 

TORQUE = 2 . 2528E+00 (N-m) , FILM POWER LOSS = 8.4929E+02 (WATT) 

AXIAL REYNOLDS NUMBER = 1.2922E+04 

CIRC. REYNOLDS NUMBERS FOR ROTOR AT SEAL ENDS = 4.2258E+03, 4.2258E+03 

DYNAMIC COEFFICIENTS ( FORCE UNIT / DISP. UNIT ) 


DISP. 

x (m) 

y (m) 

phi (RAD) 

psi (RAD) 

FORCE UNIT 

Kx 

1 . 8896E+07 

4 . 1269E+06 

- 3 . 3389E+04 

1 . 4185E+06 

N 

Ky 

-4 . 1269E+06 

1. 8896E+07 

-1.4185E+06 

-3 . 3389E+04 

N 

Kphi 

-1 . 2847E+04 

9. 8025E+04 

-3 . 8 943E+03 

1. 0298E+02 

N-m 

Kpsi 

- 9 . 8025E+04 

-1 . 2847E+04 

-1 . 02 98E+02 

-3 . 8 943E + 03 

N-m 

Bx 

2 . 1895E+04 

1. 1396E+03 

-7 . 0298E- 01 

1 . 7716E+02 

N-SEC 

By 

- 1 . 13 96E+03 

2 . 18 95E+04 

- 1 . 7716E+02 

-7 . 0298E- 01 

N-SEC 

Bphi 

- 1 . 93 73E- 01 

6 . 8162E+01 

5 . 4617E- 01 

1 . 3431E- 02 

N-m-SEC 

Bpsi 

- 6 . 8162E+01 

-1. 9373E-01 

- 1 . 3431E- 02 

5 . 4617E- 01 

N-m-SEC 

Ax 

3 . 0199E+00 

-1. 1981E-02 

2.1688E-04 

1 . 8092E-03 

N-SEC2 

Ay 

1 . 1981E- 02 

3 . 0199E+00 

- 1 . 8 092E- 03 

2 . 1688E- 04 

N-SEC2 

Aphi 

6 . 0437E-05 

4. 9833E-04 

3 . 5919E-05 

1 . 1575E- 06 

N-m-SEC2 

Apsi 

-4 . 9833E-04 

6 . 0437E-05 

-1 . 1575E- 06 

3 . 5919E- 05 

N-m- SEC2 
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( CASE 6 

) 

Childs finite 

length ; 

solution, RPM0= 

RPM/ 2 , L/D=l 


&INPUTS 








TITLE 

= 

’Childs finite 

length 

solution, RPM0 

=RPM/ 2 , L/D=l ' 

i 

I FACE 

= 

0 

ISIUN 

= 

1 



I GROT 

= 

0 

NO I 

= 

0 

I FLOW = 

0 

R0 

= 

7 . 6200E- 02 

EL 

= 

1 . 524 0E- 01 

C 

1. 9050E-04 

RPM 

= 

3 . 6000E+03 

RPM0 

= 

1 . 8000E+03 

RPMD 

0.0000E+00 

PLEG 

= 

3 . 4400E+06 

PRIG 

= 

0 . 0000E+00 

FZD 

0 . 0000E+00 

vise 

= 

1 . 2950E- 03 

DENS 

= 

1 . 0000E+03 



EMA 

= 

-2 . 5000E- 01 

ENA 

= 

7 . 9100E- 02 



EMB 

= 

-2 . 5000E- 01 

ENB 

= 

7 . 9100E- 02 



HTAP 

= 

0.0000E+00 

HBRL 

= 

0 . 0000E+00 



TOLH 

= 

1.0000E-04 

TOLV 

= 

1. 0000E-05 

DUT 

1. 0000E-06 

I HOME 

= 

0 

NITH 

= 

10 

NITV 

30 

NREG 

= 

1 

NRSUB 

= 

200 



ELFR 

= 

1 . 0000E+00 






ZET 

= 

1. 0000E-01 






ALP I 

= 

0.0000E+00 






BETI 

= 

0.0000E+00 






DELT 

= 

0 . 0000E+00 






NSG 

= 

0 






ZETG 

= 

0.0000E+00 







/ 


CYLINDRICAL SEAL 

LENGTH, DIAMETER, CLEARANCE = 1.5240E-01, 1.5240E-01, 1.9050E-04 (m) 

ROTOR, SWIRL AND DIST. SPEEDS = 3.6000E+03, 1.8000E+03, 0.0000E+00 (RPM) 

PRESSURE AT START, END AXIAL BOUNDARIES = 3.4400E+06', O.OOOOE+OO (Pa) 

VISCOSITY = 1 . 2950E-03 (Pa-SEC) , DENSITY = 1.0000E+03 (Kg/m3) 

ERROR CODE = 0, ITERATIONS IN PRIMARY FLOW = 4 

FLOW = 1 . 7711E-03 (m**3/SEC) 

TORQUE = 6 . 9241E+00 (N-m) , FILM POWER LOSS = 2.6103E+03 (WATT) 

AXIAL REYNOLDS NUMBER = 5.7131E+03 

CIRC. REYNOLDS NUMBERS FOR ROTOR AT SEAL ENDS = 4.2258E+03, 4.2258E+03 

DYNAMIC COEFFICIENTS ( FORCE UNIT / DISP. UNIT ) 


DISP. 

x (m) 

y (m) 

phi (RAD) 

psi (RAD) 

FORCE UNIT 

Kx 

1. 0794E+07 

9 . 1778E+07 

- 1 . 5134E+06 

1 . 5966E+07 

N 

Ky 

- 9 . 1778E+07 

1. 0794E+07 

- 1 . 5966E+07 

-1 . 5134E+06 

N 

Kphi 

-4.6 942E+05 

6 . 5852E+05 

-4 . 7345E+04 

5 . 3912E+04 

N-m 

Kpsi 

-6 . 5852E+05 

-4.6 942E+05 

- 5 . 3 912E+04 

-4 . 7345E+04 

N-m 

Bx 

4 . 8718E+05 

1 . 0293E+05 

1 . 1505E+02 

8 . 0173E+03 

N-SEC 

By 

-1. 0293E+05 

4 . 8 718E + 05 

-8 . 0173E+03 

1 . 1505E+02 

N-SEC 

Bphi 

-6 . 5590E+01 

2 . 4 967E+03 

2 . 8624E+02 

5 . 1787E+01 

N-m-SEC 

Bpsi 

-2 . 4967E+03 

-6 . 5 590E+01 

-5 . 1787E+01 

2 . 8624E+02 

N-m-SEC 

Ax 

2 . 7261E+02 

-2 . 1579E+00 

-8 . 7333E- 02 

-2 . 8743E- 01 

N-SEC2 

Ay 

2 . 1579E+00 

2 . 7261E+02 

2 . 8743E-01 

-8 . 7333E- 02 

N-SEC2 

Aphi 

4 . 7908E-02 

1 . 6434E- 01 

1 . 3 7 03E- 01 

- 1 . 7189E- 03 

N-m-SEC2 

Apsi 

- 1 . 6434E- 01 

4 . 7908E- 02 

1 . 7189E-03 

1 . 3703E- 01 

N-m-SEC2 
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( CASE 7 

) 

Childs finite 

length solution, RPM0=0, 

L/D= . 2 


&INPUTS 








TITLE 

= 

'Childs finite 

length 

solution, RPM0=0 

, L/D= . 2 ' 


I FACE 

= 

0 

ISIUN 

= 

1 



I GROT 

= 

0 

NO I 

= 

0 

I FLOW = 

0 

R0 

= 

7 . 6200E- 02 

EL 

= 

3 . 0480E-02 

C 

1 . 9050E- 04 

RPM 

= 

3 . 6000E+03 

RPM0 

= 

0 . 0000E+00 

RPMD 

0 . 0000E+00 

PLEG 

= 

3 . 4400E+06 

PRIG 

= 

0 . 0000E+00 

FZD 

0 . 0000E+00 

vise 

= 

1 . 2950E- 03 

DENS 

= 

1 . 0000E+03 



EMA 

= 

-2 . 5000E-01 

ENA 

= 

7 . 9100E- 02 



EMB 

= 

-2 . 5000E- 01 

ENB 

= 

7 . 9100E- 02 



HTAP 

= 

0 . 0000E+00 

HBRL 

= 

0 . 0000E+00 



TOLH 

= 

1. 0000E-04 

TOLV 

= 

1.0000E-05 

DUT 

1.0000E-06 

I HOME 

= 

0 

NITH 

= 

10 

NITV 

30 

NREG 

= 

1 

NRSUB 

= 

200 



ELFR 

= 

1. 0000E+00 






ZET 

= 

1.0000E-01 






ALP I 

= 

0 . 0000E+00 






BETI 

= 

0 . 0000E+00 






DELT 

= 

0 . 0000E+00 






NSG 

= 

0 






ZETG 

/ 

= 

0 . 0000E+00 






CYLINDRICAL 

SEAL 







LENGTH, DIAMETER, CLEARANCE = 3.0480E-02, 1.5240E-01, 1.9050E-04 (m) 

ROTOR, SWIRL AND DI ST. SPEEDS = 3.6000E+03, 0.0000E+00, 0 . 0000E+00 (RPM) 

PRESSURE AT START, END AXIAL BOUNDARIES = 3.4400E+06, 0.0000E+00 (Pa) 

VISCOSITY = 1 . 2950E-03 (Pa-SEC) , DENSITY = 1.0000E+03 (Kg/m3) 

ERROR CODE = 0, ITERATIONS IN PRIMARY FLOW = 4 

FLOW = 3 . 9890E-03 (m**3/SEC) 

TORQUE = 3 . 6677E+00 (N-m) , FILM POWER LOSS = 1.3827E+03 (WATT) 

AXIAL REYNOLDS NUMBER = 1.2867E+04 

CIRC. REYNOLDS NUMBERS FOR ROTOR AT SEAL ENDS = 8.4516E+03, 5.3357E+03 


DYNAMIC COEFFICIENTS ( FORCE UNIT / DISP . UNIT ) 


DISP. 

x (m) 

y (m) 

phi (RAD) 

psi (RAD) 

FORCE UNIT 

Kx 

1 . 8583E+07 

-3 . 0274E+05 

-3 . 0756E+04 

1 . 4106E+06 

N 

Ky 

3 . 0274E+05 

1 . 8583E+07 

- 1 . 4106E+06 

-3 . 0756E+04 

N 

Kphi 

2 . 1014E+03 

9 . 6616E+04 

-3 . 8676E+03 

-1 . 6999E+01 

N-m 

Kpsi 

- 9 . 6616E+04 

2 . 1014E+03 

1 . 6999E+01 

-3 . 8676E+03 

N-m 

Bx 

2 . 1892E+04 

8 . 5184E+02 

1 . 6580E+00 

1 . 7615E+02 

N-SEC 

By 

-8 . 5184E+02 

2 . 1892E+04 

-1 . 7615E+02 

1 . 6580E + 00 

N-SEC 

Bphi 

- 1 . 0540E+00 

6 . 8066E+01 

5 . 4479E- 01 

1 . 9870E- 02 

N-m-SEC 

Bpsi 

- 6 . 8 066E+01 

- 1 . 0540E+00 

-1 . 9870E-02 

5 . 4479E- 01 

N-m-SEC 

Ax 

3 . 0025E+00 

- 8 . 3 743E- 02 

- 5 . 4 944E- 04 

1 . 8925E- 03 

N-SEC2 

Ay 

8 . 3743E- 02 

3 . 0025E+00 

-1 . 8925E-03 

- 5 . 4944E- 04 

N-SEC2 

Aphi 

3 . 7677E- 04 

4 . 5070E- 04 

3.5698E-05 

-2 . 5857E- 06 

N-m-SEC2 

Apsi 

-4 . 5070E- 04 

3 . 7677E- 04 

2 . 5857E- 06 

3 . 5698E- 05 

N-m-SEC2 
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( CASE 8 

) 

Childs finite 

length i 

solution, RPM0=0 , 

l/d=i 


&INPUTS 








TITLE 

= 

'Childs finite 

length 

solution, RPM0=0 

, l/d=i ' 

i 

I FACE 

= 

0 

ISIUN 

= 

1 



I GROT 

= 

0 

NOI 

= 

0 

IFLOW 

= 0 

R0 

= 

7 . 6200E- 02 

EL 

= 

1 . 5240E- 01 

C 

= 1 . 9050E- 04 

RPM 

= 

3 . 6000E+03 

RPM0 

= 

0 . 0000E+00 

RPMD 

= 0 . 0000E+00 

PLEG 

= 

3 . 4400E+06 

PRIG 

= 

0 . 0000E+00 

FZD 

= 0 . 0000E+00 

vise 

= 

1 . 2950E- 03 

DENS 

= 

1. 0000E+03 



EMA 

= 

-2 . 5000E- 01 

ENA 

= 

7 . 9100E- 02 



EMB 

= 

-2 . 5000E- 01 

ENB 

= 

7 . 9100E- 02 



HTAP 

= 

0 . 0000E+00 

HBRL 

= 

0 . 0000E+00 



TOLH 

= 

1. 0000E-04 

TOLV 

= 

1. 0000E-05 

DUT 

= 1.0000E-06 

I HOME 

= 

0 

NITH 

= 

10 

NITV 

= 30 

NREG 

= 

1 

NRSUB 

= 

200 



ELFR 

= 

1 . 0000E+00 






ZET 

= 

1.0000E-01 






ALPI 

= 

0 . 0000E+00 






BETI 

= 

0 . 0000E+00 






DELT 

= 

0 . 0000E+00 






NSG 

= 

0 






ZETG 

= 

0 . 0000E+00 







/ 


CYLINDRICAL SEAL 

LENGTH, DIAMETER, CLEARANCE = 1.5240E-01, 1.5240E-01, 1.9050E-04 (m) 

ROTOR, SWIRL AND DIST. SPEEDS = 3.6000E+03, 0.0000E+00, 0.0000E+00 (RPM) 

PRESSURE AT START, END AXIAL BOUNDARIES = 3.4400E+06, 0.0000E+00 (Pa) 

VISCOSITY = 1 . 2 950E- 03 (Pa-SEC) , DENSITY = 1.0000E+03 (Kg/m3) 

ERROR CODE = 0, ITERATIONS IN PRIMARY FLOW = 4 
FLOW = 1 . 7673E-03 (m**3/SEC) 

TORQUE = 7 . 8249E+00 (N-m), FILM POWER LOSS = 2 . 9499E+03 (WATT) 

AXIAL REYNOLDS NUMBER = 5.7007E+03 

CIRC. REYNOLDS NUMBERS FOR ROTOR AT SEAL ENDS = 8.4516E+03, 4.2259E+03 

DYNAMIC COEFFICIENTS ( FORCE UNIT / DISP. UNIT ) 


DISP. 

x (m) 

y (m) 

phi (RAD) 

psi (RAD) 

FORCE UNIT 

Kx 

1 . 3250E+07 

7 . 5182E+07 

- 1 . 8882E+06 

1 . 5946E+07 

N 

Ky 

-7 . 5182E+07 

1 . 3250E+07 

- 1 . 5946E+07 

-1 . 8882E+06 

N 

Kphi 

-3 . 3283E+04 

7 . 0821E+05 

-4 . 8210E+04 

4 . 1331E+04 

N-m 

Kpsi 

-7 . 0821E+05 

-3 . 3283E+04 

-4 . 1331E+04 

-4 . 8210E+04 

N-m 

Bx 

4 . 8952E+05 

8 . 9262E+04 

- 1 . 4507E+02 

7 . 903 0E+03 

N-SEC 

By 

-8 . 9262E+04 

4 . 8952E+05 

- 7 . 903 0E+03 

- 1 . 4507E+02 

N-SEC 

Bphi 

2 . 5693E+02 

2 . 563 9E+03 

2 . 8825E+02 

4 . 4699E+01 

N-m-SEC 

Bpsi 

-2 . 5639E+03 

2 . 56 93E+02 

-4 . 46 99E+01 

2 . 8825E+02 

N-m-SEC 

Ax 

2 . 7214E+02 

-3 . 2 53 8E+00 

-1 . 5739E-01 

-2 . 566 0E- 01 

N-SEC2 

Ay 

3 . 2538E+00 

2 . 72 14E + 02 

2 . 5660E- 01 

-1.5739E-01 

N-SEC2 

Aphi 

6 . 4178E- 02 

1. 5751E-01 

1 . 3636E- 01 

-3 . 2263E-03 

N-m-SEC2 

Apsi 

-1 . 5751E-01 

6 . 4178E- 02 

3 . 2263E- 03 

1 . 3636E- 01 

N-m-SEC2 
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6. Operating Environment 


The computer code SPIRALI has been written to run under a variety of operating environments with very little 
modification. The extensions to ANSI Standard FORTRAN 77 that have been used are discussed at the start 
of Section 3. Executable versions of SPIRALI have been compiled with Version 5.1 of the Microsoft 
FORTRAN compiler and tested for use on IBM PC compatible computers with 80x87 floating point 
coprocessors. The source program listed in Appendix A has been compiled with the command: 

FL /AH /Gtl 024 /c SPIRALI.FOR 

and linked with the appropriate large model, floating point libraries provided by Microsoft to run under either 
DOS or OS2 operating systems. 

The memory required to run the program is controlled by the values of the parameters NDZ and NDREG in 
the one line parameter statement that appears in the main program and many of the supporting sub-programs. 
It currently has the form 


PARAMETER (NDZ=201,NDREG=21) 

permitting values of NRSUB of up to 200 in each of 21 regions and using approximately 500k of available 
memory. 

The DOS version of the program should run on any IBM PC or compatible with a math coprocessor, sufficient 
memory and DOS 3.0 or higher. Mass storage requirements are minimal (a 360k diskette drive should be 
sufficient). The OS2 version will require sufficient additional resources to run OS2 (version 1 .2 or higher). 
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7. Error Messages 


1 Initial velocity computation diverged. 

Occurs when the Newton iterative procedure for finding an initial scale factor for the velocity diverges. 
Check to see that reasonable values are used for C, DUT and NITV. 

2 Primary flow computation diverged. 

Occurs when the Newton iterative procedure for finding the transverse velocity at the inlet diverges 
in the primary flow solution. This error is likely to occur when transverse inertia effects are small but 
not neglected and the grid is too course. Try using a finer grid, particularly near the inlet. If that fails 
set NOI=1 . Check to see that reasonable values are used for DUT and NITV. 

3 Matrix inversion error encountered in second order solution. 

This error occurs when the derivative matrix used in the semi-implicit solution algorithm becomes 
singular. If inputs are reasonable this error can be corrected by using a finer grid or setting NOI=1. 

4 Spiral groove local flow computation diverged. 

Occurs when the Newton-Raphson procedure for determining the global turbulent shear functions in 
grooved regions diverges. Check to see that the input spiral groove parameters are correct. If so, 
using a finer grid should eliminate this problem. 

5 Face seal axial load iteration diverged. 

This error can readily occur when IHOME=1 and a bad initial guess is used for the starting value of 
C. Try running with IHOME=0 to find a good starting value of C. Also make sure that the grid is fine 
enough to avoid excessive truncation error and that you have used a sufficiently large value for NITH. 

6 Negative stiffness or film thickness in axial load iteration. 

Causes and remedies same as error 5 above except changing NITH won't help. 

7 Wrong inlet boundary was used with transverse inertia included. 

Change value of IFLOW or set NOI=1 . 

8 Illegal length, clearance, viscosity, pressure or speed encountered. 

Check inputs and correct error. 

9 Maximum number of allowable grid points exceeded. 

Reduce values of NRSUB. 
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10 Maximum number of allowable regions exceeded. 

Reduce NREG. 

11 Sum of length fractions are not equal to 1. 

Adjust values of ELFR so that they add up to 1 . 
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APPENDIX A 


Application to Parallel and Helical Grooved Seals 
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Parallel groove geometries can be analyzed by using the separated region option of SPIRALI. Analysis 
was conducted of the parallel groove and stator 1 of Reference 1 by treating each groove and ridge as a 
separate region. The clearance and groove geometry are shown on Figure 1. 

A total of 19 regions were used consisting of nine grooves, eight lands between grooves, an inlet land and 
an exit land. (The parameter statements in SPIRALI governing the maximum number of regions were 
enlarged appropriately). The inlet and exit land lengths were enlarged slightly to make the sum of the 
lengths of the individual regions equal to the total seal length (50.8 mm). The seal radius is 50.8 mm. The 
viscosity and density are 1.54 x 104 Pa-s and 1570 kg/m 3 , respectively. The program default values for 
flow parameters of m = -0.25, n = 0.0791 (Hirs coefficients for smooth surfaces as given by Blasius theory), 
and loss coefficient of £ = 0 were used. Although it was stated in Reference 1 that no intentional swirl was 
provided, a swirl speed of 25% of the rotor speed was used in all computations to characterize the swirl 
that might naturally be present. 

Comparison of Flow Rates and Force Coefficients 

Experimental flow rates are presented in terms of a discharge coefficient, C D , defined by Equation (3) of 
Reference 1 . The value of C D ' 1/2 for stator 1 is shown in Figure 6 of Reference 1 to have a constant value 
of approximately 0.5 for all measured pressures. A comparison between that value and those predicted 
by SPIRALI is shown in Figure 2. The agreement is good and the results are predicted to be relatively 
insensitive to the rotating speed over the range of conditions covered, as indicated in Reference 1. 

The dynamic coefficients presented in Table 1 of Reference 1 are computed from measured force 
coefficients based on the assumption that the coefficients are independent of speed. For purposes of 
comparison, the measured force coefficients, f B and f r , were retrieved with the use of Equation 6 in 
Reference 1 . The seal pressure drops, AP, were calculated from the axial Reynolds number, R a , using 
Equation (3) of Reference 1 with C D ' 1/2 = 0.5. Values of f e versus AP obtained in this manner and those 
predicted by SPIRALI appear to be in good agreement, as shown in Figure 3. The predicted variation in 
f r with rotating speed (Figure 4) is somewhat weaker than that obtained by experiment. Discrepancies in 
f r can, in part, be due to sensitivity to the input values to SPIRALI provided for the friction factor and loss 
coefficients. Also, the effective mass, M ef , in Table 1 of Reference 1 shows a significant jump at a 
Reynolds number of 200,000 that may be suspect. 
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Effective stiffness values are compared in Figure 5. The experimental stiffnesses are based on 
extrapolations to zero speed. Stiffness is assumed to be caused principally by the pressure gradient and 
inlet Bernoulli effect which causes a circumferential pressure gradient and positive stiffness if the rotor goes 
eccentric. As Figure 5 indicates, the predicted values approach the experimental values as the speed is 
reduced. The predicted decrease in K with increasing rotating speed is to be expected due to the Bernoulli 
effect associated with fluid rotation. A corresponding comparison of C ef with C - k/oo is given in Figure 6. 
In both cases, the agreement improves as the rotational speed is decreased and the physical interpretation 
of the data becomes valid. 

Effect of Local Pressure Discontinuities on Predicted Flow 

As described above, the effects of pressure jumps at the grooves are accounted for by SPIRALI in the prediction 
of flow and force coefficients. An estimate of these terms can be obtained by comparing the flow rate for circular 
grooves with that which would be predicted to occur if the effects of pressure jumps at the groove interfaces 
were neglected. This comparison is made at 1000 rpm in Figure 7. The large discrepancies between the two 
curves indicate a significant effect of local inertia that must be accounted for. 

Check cases were also run against the helical groove results of Reference 1 . Figure 8 shows comparisons 
of the flow coefficient for varying groove angles at a loss coefficient equal to 1 , and the theory predicts 
results accurately. Cross-coupled force comparisons are shown on Figure 9. Cross-coupled effects are 
important because they are a measure of the stability of the seal and in some cases cross-coupled forces 
cause subsynchronous whirling of the rotor. The comparative results are very good. It is noted that cross- 
coupled forces were compared rather than the stiffnesses indicated in Reference 1 . The experimental 
procedure measures forces and computes stiffness and damping based on the assumption that the cross- 
coupled coefficients vary directly with speed. The assumption of zero cross coupling at zero speed is 
suspect for helical grooves because the grooves apply a pressure-driven tangential flow at zero speed 
opposite the direction of rotation. As shown on Figure 10, the comparisons of direct stiffness, K ef , were not 
as good as the parameters presented above. This may be due to the experimental assumptions or code 
assumptions. 

The net results of these studies is that SPIRALI is an effective tool for analyzing many parallel and helical 
groove seals. The principal limitation is the bulk flow model that does not treat flow variations across gaps. 
If large gaps are used, then more sophisticated CFD analysis is required. 
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Figure 1 Parallel-groove pressure breakdown seal 
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Figure 2 Parallel-groove seal flow coefficient 
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Figure 3 Tangential force coefficients 
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Figure 4 Normal force coefficients 
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Comparison between K and K ef at Various Rotating Speeds 
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Figure 5 Parallel-groove seal effective stiffness 
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Figure 6 Comparison between C - k/co and C ef at various rotating speeds 
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Figure 7 Effect of local pressure discontinuities on predicted axial flow rates 
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Figure 8 Flow coefficient: SPIRALI compared to reference 1 
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